Astrophysical Accretion as an Analogue Gravity Phenomena 
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Inspite of the remarkable resemblance in between a black hole and an ordinary thermodynamic 
system, black holes never radiate according to the classical laws of physics. The introduction of 
quantum effects radically changes the scenario. Black holes radiate due to quantum effects. Such 
radiation is known as Hawking radiation and the corresponding radiation temperature is referred 
as the Hawking temperature. Observational manifestation of Hawking effect for astrophysical black 
holes is beyond the scope of present day's experimental techniques. Also, Hawking quanta may 
posses trans-Planckian frequencies, and physics beyond the Planck scale is not well understood. The 
above mentioned difficulties with Hawking effect were the motivations to search for an analogous 
version of Hawking radiation, and the theory of acoustic/analogue black holes were thus introduced. 

Classical black hole analogues (alternatively, the analogue systems) are fluid dynamical analogue 
of general relativistic black holes. Such analogue effects may be observed when the acoustic pertur- 
bation (sound waves) propagates through a classical dissipation-less transonic fluid. The acoustic 
horizon, which resembles the actual black hole event horizon in many ways, may be generated at 
the transonic point in the fluid flow. Acoustic horizon emits acoustic radiation with quasi thermal 
phonon spectra, which is analogous to the actual Hawking radiation. 

Transonic accretion onto astrophysical black holes is a very interesting example of classical ana- 
logue system found naturally in the Universe. An accreting black hole system as a classical analogue 
is unique in the sense that only for such a system, both kind of horizons, the electromagnetic and the 
acoustic (generated due to transonicity of accreting fluid) are simultaneously present in the same 
system. Hence accreting astrophysical black holes are the most ideal candidate to study theoreti- 
cally and to compare the properties of these two different kind of horizons. Such a system is also 
unique in the aspect that accretion onto the black holes represents the only classical analogue system 
found in the nature so far, where the analogue Hawking temperature may exceed the actual Hawking 
temperature. In this review article, it will be demonstrated that, in general, the transonic accretion 
in astrophysics can be considered as an example of the classical analogue gravity model. 
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I. BLACK HOLES 

Black holes are the vacuum solutions of Einstein's field equations in general relativity. Classically, a black hole is 
conceived as a singularity in space time, censored from the rest of the Universe by a mathematically defined one way 
surface, the event horizon. Black holes are completely characterized only by three externally observable parameters, 
the mass of the black hole Mbh, the rotation (spin) Jbh and charge Qbh- All other informations about the matter 
which formed the black hole or is falling into it, disappear behind the event horizon, are therefore permanently 
inaccessible to the external observer. Thus the space time metric defining the vacuum exterior of a classical black 
hole is characterized by Mbh, Jbh and Qbh only. The most general family of black hole solutions have non zero 
values of Mbh, Jbh and Qbh (rotating charged black holes), and are known as the Kerr-Newman black holes. The 
following table classifies various categories of black hole solutions according to the value of Mbh, Jbh and Qbh- 



Types of Black Hole 


Mass 


Angular Momentum 


Charge 


Kerr-Newman 
(Newman et. al. 1965) 


Mbh > 


Jbht^O 


Qbht^O 


Kerr (Kerr 1963) 


Mbh > 


Jbh^O 


Qbh=0 


Reissner-N ordstrom 

(Reissner 1916; Weyl 1917; Nordstrom 1918) 


Mbh > 


Jbh=0 


Qbht^O 


Schwarzs child 
(Schwarzschild 1916) 


Mbh > 


Jbh=0 


Qbh=0 
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Table 1: Classification of black holes according to the value of its Mass, angular momentum and charge. 

The Israle-Carter- Robinson theorem (Israle 1967; Carter 1971; Robinson 1975), when coupled with Price's conjec- 
ture (Price 1972), ensures that any object with event horizon must rapidly settles down to the Kerr metric, radiating 
away all its irregularities and distortions which may deviate them from the black hole solutions exactly described by 
the Kerr metric. 

In astrophysics, black holes are the end point of gravitational collapse of massive celestial objects. The Kerr- 
Newman and the Reissner-Nordstrom black hole solutions usually do not play any significant role in astrophysical 
context. Typical astrophysical black holes are supposed to be immersed in an charged plasma environment. Any 
net charge Qbh will thus rapidly be neutrilized by teh ambient magnetic field. The time scale for such charge 
relaxation would be roughly of the order of (Mbh/Mq) fisec [Mq being the mass of the Sun, see, e.g., Hughes 2005 
for further details), which is obviously far shorter compared to the rather long timescalc relevant to observing most 
of the properties of the astrophysical black holes. Hence the Kerr solution provides the complete description of most 
stable astrophysical black holes. However, the study of Schwarzschild black holes, although less general compared to 
the Kerr type holes, is still greatly relevant in astrophysics. 

Astrophysical black holes may be broadly classified into two categories, the stellar mass [Mbh^ a few Mq), and 
super massive {Mbh^^^^^q) black holes. While the birth history of the stellar mass black holes is theoretically 
known with almost absolute certainty (they are the endpoint of the gravitational collapse of massive stars), the 
formation scenario of the supermassive black hole is not unanimously understood. A super massive black hole may 
form through the monolithic collapse of early proto-spheroid gaseous mass originated at the time of galaxy formation. 
Or a number of stellar/intermediate mass black holes may merge to form it. Also the runaway growth of a seed black 
hole by accretion in a specially favoured high-density environment may lead to the formation of super massive black 
holes. However, it is yet to be well understood exactly which of the above mentioned processes routes toward the 
formation of super massive black holes; see, e.g., Rees 1984, 2002; Haiman & Quataert 2004; and Volonteri 2006, for 
comprehensive review on the formation and evolution of super massive black holes. 

Both kind of astrophysical black holes, the stellar mass and super massive black holes, however, accrete matter 
from the surroundings. Depending on the intrinsic angular momentum content of accreting material, either spherically 
symmetric (zero angular momentum flow of matter), or axisymmetric (matter flow with non-zero finite angular 
momentum) flow geometry is invoked to study an accreting black hole system (see the excellent monographs by Frank, 
King & Raine 1992, and Kato, Fukue & Mineshige 1998, for details about the astrophysical accretion processes). We 
will get back to the accretion process in greater detail in subsequent sections. 

II. BLACK HOLE THERMODYNAMICS 

Within the framework of purely classical physics, black holes in any diffeomorphism covariant theory of gravity (where 
the field equations directly follow from the diffeomorphism covariant Lagrangian) and in general relativity, mathe- 
matically resembles some aspects of classical thermodynamic systems (Wald 1984, 1994, 2001; Keifer 1998; Brown 
1995, and references therein). In early seventies, a series of influential works (Bekenstein 1972, 1972a, 1973, 1975; 
Israel 1976; Bardeen, Carter & Hawking 1973, see also Bekenstein 1980 for a review) revealed the idea that classical 
black holes in general relativity, obey certain laws which bear remarkable analogy to the ordinary laws of classical 
thermodynamics. Such analogy between black hole mechanics and ordinary thermodynamics ( 'The Generalized Sec- 
ond Law', as it is customarily called) leads to the idea of the 'surface gravity' of black hole. [223j k, which can be 
obtained by computing the norm of the gradient of the norms of the Killing fields evaluated at the stationary black 
hole horizon, and is found to be constant on the horizon (analogous to the constancy of temperature T on a body 
in thermal equilibrium - the 'Zeroth Law' of classical thermodynamics). Also, k = can not be accomplished by 
performing finite number of operations (analogous to the 'weak version' of the third law of classical thermodynamics 
where temperature of a system cannot be made to reach at absolute zero, see discussions in Keifer 1998). It was found 
by analogy via black hole uniqueness theorem (see, e.g., Heusler 1996, and references therein) that the role of entropy 
in classical thermodynamic system is played by a constant multiple of the surface area of a classical black hole. 

III. HAWKING RADIATION 

The resemblance between the laws of ordinary thermodynamics to those of black hole mechanics were, however, 
initially regarded as purely formal. This is because, the physical temperature of a black hole is absolute zero (see, 
e.g. Wald 2001). Hence physical relationship between the surface gravity of the black hole and the temperature of 
a classical thermodynamic system can not be conceived. This further indicates that a classical black hole can never 
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radiate. However, introduction of quantum effects might bring a radical change to the situation. In an epoch making 
paper published in 1975, Hawking (Hawking 1975) used quantum field theoretic calculation on curved spacetime to 
show that the physical temperature and entropy of black hole does have finite non-zero value (see Page 2004 and 
Padmanabhan 2005 for intelligible reviews of black hole thermodynamics and Hawking radiation) . A classical space 
time describing gravitational collapse leading to the formation of a Schwarzschild black hole was assumed to be the 
dynamical back ground, and a linear quantum field, initially in it's vacuum state prior to the collapse, was considered 
to propagate against this background. The vacuum expectation value of the energy momentum tensor of this field 
turned out to be negative near the horizon. This phenomenon leads to the flux of negative energy into the hole. Such 
negative energy flux would decrease the mass of the bl ack h ole and would lead to the fact that the quantum state of 



the outgoing mode of the field would contain particles. [224j The expected number of such particles would correspond 



to radiation from a perfect black body of finite size. Hence the spectrum of such radiation is thermal in nature, and 
the temperature of such radiation, the Hawking temperature Th from a Schwarzschild black hole, can be computed 
as 

^" " STTkeGMBH 

where G is the universal gravitational constant, c, h and ks are the velocity of light in vacuum, the Dirac's constant 
and the Boltzmann's constant, respectively. 

The semi classical description for Hawking radiation treats the gravitational field classically and the quantized 
radiation field satisfies the d'Alembert equation. At any time, black hole evaporation is an adiabatic process if the 
residual mass of the hole at that time remains larger than the Planck mass. 



IV. TOWARD AN ANALOGY OF HAWKING EFFECT 



Substituting the values of the fundamental constants in Eq. (1), one can rewrite Th for a Schwarzschild black hole 
as (Heifer 2003): 

- 6. 2 xlO^M-^] Degree Kelvin (2) 

V ^^BH J 

It is evident from the above equation that for one solar mass black hole, the value of the Hawking temperature would 
be too small to be experimentally detected. A rough estimate shows that Th for stellar mass black holes would be 
around lO'' times colder than the cosmic microwave background radiation. The situation for super massive black hole 
will be much more worse, as ThozI/ Mbh ■ Hence Th would be a measurable quantity only for primordial black holes 
with very small size and mass, if such black holes really exist, and if instruments can be fabricated to detect them. 
The lower bound of mass for such black holes may be estimated analytically. The time-scale T (in years) over which 
the mass of the black hole changes significantly due to the Hawking's process may be obtained as (Heifer 2003): 

As the above time scale is a measure of the lifetime of the hole itself, the lower bound for a primordial hole may be 
obtained by setting T equal to the present age of the Universe. Hence the lower bound for the mass of the primordial 
black holes comes out to be around 10^^ gm. The size of such a black hole would be of the order of 10"^'^ cm and the 
corresponding Th would be about lO^^K, which is comparable with the macroscopic fluid temperature of the freely 
falling matter (spherically symmetric accretion) onto an one solar mass isolated Schwarzschild black hole (see section 
12.1 for further details). However, present day instrumental technique is far from efficient to detect these primordial 
black holes with such an extremely small dimension, if such holes exist at all in first place. Hence, the observational 
manifestation of Hawking radiation seems to be practically impossible. 

On the other hand, due to the infinite redshift caused by the event horizon, the initial configuration of the 
emergent Hawking Quanta is supposed to possess trans-Planckian frequencies and the corresponding wave lengths are 
beyond the Planck scale. Hence, low energy effective theories cannot self consistently deal with the Hawking radiation 
(see, e.g., Parentani 2002 for further details). Also, the nature of the fundamental degrees of freedom and the physics 
of such ultra short distance is yet to be well understood. Hence, some of the fundamental issues like the statistical 
meaning of the black hole entropy, or the exact physical origin of the out going mode of the quantum field, remains 
unresolved (Wald 2001). 

Perhaps the above mentioned difhculties associated with the theory of Hawking radiation served as the principal 
motivation to launch a theory, analogous to the Hawking's one, effects of which would be possible to comprehend 
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through relatively more perceivable physical systems. The theory of analogue Hawking radiation opens up the possi- 
bility to experimentally verify some basic features of black hole physics by creating the sonic horizons in the laboratory. 
A nu mber of works have been carried out to formulate the condensed matter or optical analogue of event horizons 
[225*1 . The theory of analogue Hawking radiation may find important uses in the fields of investigation of quasi- normal 
modes (Berti, Cardoso & Lemos 2004; Cardoso, Lemos & Yoshida 2004), acoustic super-radiance (Basak & Majum- 
dar 2003; Basak 2005; Lepe & Saavedra 2005; Slatyer, & Savage 2005; Cherubini, Federici & Succi 2005; Kim, Son, 
& Yoon 2005; Choy, Kruk, Carrington, Fugleberg, Zahn, Kobes, Kunstatter & Pickering 2005; Federici, Cherubini, 
Succi & Tosi 2005), FRW cosmology (Barcelo, Liberati & Visser 2003) inflationary models, quantum gravity and 
sub-Planckian models of string theory (Parentani 2002). 

For space limitation, in this article, we will, however, mainly describe the formalism behind the classical analogue 
systems. By 'classical analogue systems' we refer to the examples where the analogue effects are studied in classical 
systems (fluids), and not in quantum fluids. In the following sections, we discuss the basic features of a classical 
analogue system. 



V. ANALOGUE GRAVITY MODEL AND THE BLACK HOLE ANALOGUE 



In recent years, strong analogies have been established between the physics of acoustic perturbations in an inhomo- 
geneous dynamical fluid system, and some kinematic features of space-time in general relativity. An effective metric, 
referred to as the 'acoustic metric', which describes the geometry of the manifold in which acoustic perturbations 
propagate, can be constructed. This effective geometry can capture the properties of curved space-time in general 
relativity. Physical models constructed utilizing such analogies are called 'analogue gravity models' (for details on 
analogue gravity models, see, e.g. the review articles by Barcelo, Liberati & Visser (2005) and Cardoso (2005), and 
the monograph by Novello, Visser & Volovik (2002)). 

One of the most significant effects of analogue gravity is the 'classical black hole analogue'. Classical black 
hole analogue effects may be observed when acoustic perturbations (sound waves) propagate through a classical, 
dissipation-less, inhomogeneous transonic fluid. Any acoustic perturbation, dragged by a supersonically moving fluid, 
can never escape upstream by penetrating the 'sonic surface'. Such a sonic surface is a collection of transonic points in 
space-time, and can act as a 'trapping' surface for outgoing phonons. Hence, the sonic surface is actually an acoustic 
horizon, which resembles a black hole event horizon in many ways and is generated at the transonic point in the fluid 
flow. The acoustic horizon is essentially a null hyper surface, generators of which are the acoustic null geodesies, i.e. 
the phonons. The acoustic horizon emits acoustic radiation with quasi thermal phonon spectra, which is analogous 
to the actual Hawking radiation. The temperature of the radiation emitted from the acoustic horizon is referred to 
as the analogue Hawking temperature. 

Hereafter, we shall use Tah to denote the analogue Hawking temperature, and Th to denote the the actual 
Hawking temperature as defined in ([T]). We shall also use the words 'analogue', 'acoustic' and 'sonic' synonymously in 
describing the horizons or black holes. Also the phrases 'analogue (acoustic) Hawking radiation/effect/temperature' 
should be taken as identical in meaning with the phrase 'analogue (acoustic) radiation/effect/temperature'. A system 
manifesting the effects of analogue radiation, will be termed as analogue system. 

In a pioneering work, Unruh (1981) showed that a classical system, relatively more clearly perceivable than a 
quantum black hole system, does exist, which resembles the black hole as far as the quantum thermal radiation is 
concerned. The behaviour of a linear quantum field in a classical gravitational field was simulated by the propagation 
of acoustic disturbance in a convergent fluid flow. In such a system, it is possible to study the effect of the reaction of 
the quantum field on it's own mode of propagation and to contemplate the experimental investigation of the thermal 
emission mechanism. Considering the equation of motion for a transonic barotropic irrotational fluid, Unruh (1981) 
showed that the scaler field representing the acoustic perturbation (i.e, the propagation of sound wave) satisfies a 
differential equation which is analogous to the equation of a massless scaler field propagating in a metric. Such a 
metric closely resembles the Schwarzschild metric near the horizon. Thus acoustic propagation through a supersonic 
fluid forms an analogue of event horizon, as the 'acoustic horizon' at the transonic point. The behaviour of the normal 
modes near the acoustic horizon indicates that the acoustic wave with a quasi-thermal spectrum will be emitted from 
the acoustic horizon and the temperature of such acoustic emission may be calculated as (Unruh 1981): 



h 

Tah ~ 



AnkB 



(4) 



Where rh represents the location of the acoustic horizon, Cs is the sound speed, u± is the component of the dynamical 
flow velocity normal to the acoustic horizon, and d/drj represents derivative in the direction normal to the acoustic 
horizon. 
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Equation ([4]) has clear resemblance with ^ and hence Tah is designated as analogue Hawking temperature and 
such quasi-thermal radiation from acoustic (analogue) black hole is known as the analogue Hawking radiation. Note 
that the sound speed Cs in Unruh's original treatment (the above equation) was assumed to be constant in space, i.e., 
an isothermal equation of state had been invoked to describe the fluid. 

Unruh's work was followed by other important papers (Jacobson 1991, 1999; Unruh 1995; Visser 1998; Bilic 
1999) A more general treatment of the classical analogue radiation for Newtonian fluid was discussed by Visser 
(1998), who considered a general barotropic, inviscid fluid. The acoustic metric for a point sink was shown to be 
conformally related to the Painleve-GuUstrand-Lemaitre form of the Schwarzschild metric (Painleve 1921; GuUstrand 
1922; Lemaitre 1933) and a more general expression for analogue temperature was obtained, where unlike Unruh's 
original expression the speed of sound was allowed to depend on space coordinates. 

In the analogue gravity systems discussed above, the fluid flow is non-relativistic in flat Minkowski space, whereas 
the sound wave propagating through the non-relativistic fluid is coupled to a curved pseudo-Riemannian metric. This 
approach has been extended to relativistic fluids (Bilic 1999) by incorporating the general relativistic fluid dynamics. 

In subsequent sections, we will pedagogically develop the concept of the acoustic geometry and related quantities, 
like the acoustic surface gravity and the acoustic Hawking temperature. 



VI. CURVED ACOUSTIC GEOMETRY IN A FLAT SPACE-TIME 



Let denote the velocity potential describing the fluid flow in Newtonian space-time, i.e. let u — —Vi/j, where u is 
the velocity vector describing the dynamics of a Newtonian fluid. The speciflc enthalpy ft, of a barotropic Newtonian 
fluid satisfles V/i = (l/p)Vp, where p and p are the density and the pressure of the fluid. One then writes the Euler 
equation as 

-5tV + /i+^(VV')' + $ = 0, (5) 

where $ represents the potential associated with any external driving force. Assuming small fluctuations around some 
steady background po,Po and V'o, one can linearize the continuity and the Euler equations and obtain a wave equation 
(see Landau & Lifshitz 1959, and Visser 1998, for further detail). 
The continuity and Euler's equations may be expressed as: 

^+V-ipu) (6) 



du 



du 
'dt 



-Wp + F 



(7) 



with F being the sum of all external forces acting on the fluid which may be expressed in terms of a potential 

F = -/9V$, (8) 

Euler's equation may now be recast in the form 



^ux (V X u) Vp - V -W^ + $ 

ot p \2 



(9) 



Next we assume the fluid to be inviscid, irrotational, and barotropic. Introducing the specific enthalpy h, such that 

P 

and the velocity potential ip for which u — — Vt/^, Eq. ^ may be written as 



(10) 



^ + ft+i(VV.)2 + $ = 



(11) 



One now linearizes the continuity and Euler's equation around some unperturbed background flow variables po, 
Po, ^0- Introducing 



p = P0 + epi+O{e^), p = po + epi+0{e'^), 



(12) 
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from the continuity equation we obtain 



dt 



+ V • (po uo) = 0; 



dpi 
dt 



+ V • (pi Mo + Po ui) = 0. 



Equation (fTU]) implies 



Pi 



dh pi 
dp Po 



Using this the Unearized Euler equation reads 



dipt) 



1 



dtl-Jo , Pi 



^ + ho + 7t(VVo)' + $ = 0; + i^-uo- VVi = 0. 



at 



9i pi 



Re- arrangement of the last equation together with the barotropic assumption yields 

dp dp ,„ , ^ „ , > 
Pi = ^ Pi = ^ Po [dtVi + Uo ■ yipi) ■ 

Substitution of this into the linearized continuity equation gives the sound wave equation 



d_ 

dt 



dp (diJi 

g-^Po (^^+^o-VV.i 



dp 



Po Vt/ii - Po Wo -wr + • VV' 



dp 



dt 



Next, we define the local speed of sound by 



dp/ dp, 



where the partial derivative is taken at constant specific entropy. With help of the 4x4 matrix 

-/ 



where / is the 3x3 identity matrix, one can put Eq. ^T7\ to the form 

d^ ir-d^ij') = 0. 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



Equation (|^ describes the propagation of the linearized scalar potential tp-^ . The function represents the low 
amplitude fluctuations around the steady background {po,po,ipo) and thus describes the propagation of acoustic 
perturbation, .i.e. the propagation of sound waves. 

The form of Eq. ((20)) suggests that it may be regarded as a d'Alembert equation in curved space-time geometry. 
In any pseudo-Riemannian manifold the d'Alembertian operator can be expressed as (Misner, Thorne & Wheeler 
1973) 



□ 



1 



d. 



where is the determinant and g'^'^ is the inverse of the metric g^^,. Next, if one identifies 

one can recast the acoustic wave equation in the form (Visser 1998) 

1 



^-IG^^IC^"^ d,^' = 0, 



■\/— \GtJ.iy\ 

where G^i, is the acoustic metric tensor for the Newtonian fluid. The explicit form of G^^ is obtained as 



(21) 



(22) 



(23) 



G^i, = Po 



<cl-u^) 



(24) 
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The Lorentzian metric described by (j24p has an associated non-zero acoustic Riemann tensor for non-homogeneous, 
flowing fluids. 

Thus, the propagation of acoustic perturbation, or the sound wave, embedded in a barotropic, irrotational, non- 
dissipative Newtonian fluid flow may be described by a scalar d'Alembert equation in a curved acoustic geometry. 
The corresponding acoustic metric tensor is a matrix that depends on dynamical and thermodynamic variables 
parameterizing the fluid flow. 

For analogue systems discussed above, the fluid particles are coupled to the flat metric of Mankowski's space 
(because the governing equation for fluid dynamics in the above treatment is completely Newtonian), whereas the 
sound wave propagating through the non-relativistic fluid is coupled to the cttrwed pseudo-Riemannian metric. Phonons 
(quanta of acoustic perturbations) are the null geodesies, which generate the null surface, i.e., the acoustic horizon. 
Introduction of viscosity may destroy the Lorentzian invariance and hence the acoustic analogue is best observed in a 
vorticity free completely dissipation-less fluid (Visser 1998, and references therein). That is why, the Fermi superfluids 
and the Bose-Einstein condensates are ideal to simulate the analogue effects. 

The most important issue emerging out of the above discussions is that (see Visser 1998 and Barcelo, Liberati & 
Visser 2005 for further details): Even if the governing equation for fluid flow is completely non-relativistic (Newtonian), 
the acoustic fluctuations embedded into it are described by a curved pseudo-Riemannian geometry. This information 
is useful to portray the immense importance of the study of the acoustic black holes, i.e. the black hole analogue, or 
simply, the analogue systems. 

The acoustic metric (j24p in many aspects resembles a black hole type geometry in general relativity. For example, 
the notions such as 'ergo region' and 'horizon' may be introduced in full analogy with those of general relativistic 
black holes. For a stationary flow, the time translation Killing vector ^ = d/dt leads to the concept of acoustic ergo 
sphere as a surface at which G/j^C^^'' changes its sign. The acoustic ergo sphere is the envelop of the acoustic ergo 
region where ^'^ is space-like with respect to the acoustic metric. Through the equation G^u£,^£,^ = gtt = 1*^ — 0^, it is 
obvious that inside the ergo region the fluid is supersonic. The 'acoustic horizon' can be defined as the boundary of 
a region from which acoustic null geodesies or phonons, cannot escape. Alternatively, the acoustic horizon is defined 
as a time like hypersurface defined by the equation 



c?-it^ =0, 



(25) 



where u±^ is the component of the fluid velocity perpendicular to the acoustic horizon. Hence, any steady supersonic 
flow described in a stationary geometry by a time independent velocity vector field forms an ergo-region, inside which 
the acoustic horizon is generated at those points where the normal component of the fluid velocity is equal to the 
speed of sound. 

In analogy to general relativity, one also defines the surface gravity and the corresponding Hawking temperature 
associated with the acoustic horizon. The acoustic surface gravity may be obtained (Wald 1984) by computing the 
gradient of the norm of the Killing field which becomes null vector field at the acoustic horizon. The acoustic surface 
gravity k for a Newtonian fiuid is then given by (Visser 1998) 



1 

2^ 



d_ 

drj 



(26) 



The corresponding Hawking temperature is then defined as usual: 



Ta 



H 



2ttkb 



(27) 



VII. CURVED ACOUSTIC GEOMETRY IN A CURVED SPACE-TIME 

The above formalism may be extended to relativistic fluids in curved space-time background (Bilic 1999). The 
propagation of acoustic disturbance in a perfect relativistic inviscid irrotational fluid is also described by the wave 
equation of the form ([23]) in which the acoustic metric tensor and its inverse are defined as (Bilic 1999; Abraham, 
Bflic & Das 2006; Das, Bihc & Dasgupta 2006) 



G^. = £j [g^. + (1 - cDv^Vu] ; G'^- = ^ 



(i-4k 



(28) 



where p and h are, respectively, the rest-mass density and the specific enthalpy of the relativistic fluid, is the 
four- velocity, and g^j^ the background space-time metric. A (— , -|- -|- -I-) signature has been used to derive (|28p . The 
ergo region is again defined as the region where the stationary Killing vector ^ becomes spacelike and the acoustic 
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horizon as a timelike hypersurface the wave velocity of which equals the speed of sound at every point. The defining 
equation for the acoustic horizon is again of the form (|25p in which the three-velocity component perpendicular to 
the horizon is given by 

= , , (29) 

where ry'^ is the unit normal to the horizon. For further details about the propagation of the acoustic perturbation, 
see Abraham, Bilic & Das 2006. 

It may be shown that, the discriminant of the acoustic metric for an axisymmetric flow 

^ — ^t4> ~ GuG^^, (30) 

vanishes at the acoustic horizon. A supersonic flow is characterized by the condition 2? > 0, whereas for a subsonic 
flow, V < (Abraham, Bilic & Das 2006). According to the classiflcation of Bercelo, Liberati, Sonego & Visser (2004), 
a transition from a subsonic (T> < 0) to a supersonic (T> > 0) flow is an acoustic black hole, whereas a transition from 
a supersonic to a subsonic flow is an acoustic white hole. 

For a stationary configuration, the surface gravity can be computed in terms of the Killing vector 

- e + ^^"^ (31) 

that is null at the acoustic horizon. Following the standard procedure (Wald 1984; Bilic 1999) one finds that the 
expression 

«X^-^G^'^r;,|-(G„^x"x/?) (32) 

holds at the acoustic horizon, where the constant k is the surface gravity. From this expression one deduces the 
magnitude of the surface gravity as (see Bilic 1999; Abraham, Bilic & Das 2006; Das, Bilic & Dasgupta 2006 for 
further details) 



(W - Cs) 



(33) 



VIII. QUANTIZATION OF PHONONS AND THE HAWKING EFFECT 

The purpose of this section (has been adopted from Das, Bilic & Dasgupta 2006) is to demonstrate how the 
quantization of phonons in the presence of the acoustic horizon yields acoustic Hawking radiation. The acoustic 
perturbations considered here are classical sound waves or phonons that satisfy the massless wave equation in curved 
background, i.e. the general relativistic analogue of ((23)) . with the metric G^i, given by ((28l) . Irrespective of the 
underlying microscopic structure, acoustic perturbations are quantized. A precise quantization scheme for an analogue 
gravity system may be rather involved (Unruh & Schiitzhold 2003). However, at the scales larger than the atomic 
scales below which a perfect fluid description breaks down, the atomic substructure may be neglected and the field 
may be considered elementary. Hence, the quantization proceeds in the same way as in the case of a scalar field in 
curved space (Birrell & Davies 1982) with a suitable UV cutoff for the scales below a typical atomic size of a few A. 

For our purpose, the most convenient quantization prescription is the Euclidean path integral formulation. Con- 
sider a 2-|-l-dimensional axisymmetric geometry describing the fluid flow (since we are going to apply this on the 
equatorial plane of the axisymmetric black hole accretion disc, see section 13 for further details). The equation of 
motion (|23p with follows from the variational principle applied to the action functional 

S[ip]= J dtdrd(f>y/^G^"'df,ipd,.ip. (34) 

We define the functional integral 

Z = J 2?(/pe-^^t^l, (35) 

where S'e is the Euclidean action obtained from by setting t = ir and continuing the Euclidean time r from 
imaginary to real values. For a field theory at zero temperature, the integral over t extends up to infinity. Here, 
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owing to the presence of the acoustic horizon, the integral over r will be cut at the inverse Hawking temperature 27r/K 
where n denotes the analogue surface gravity. To illustrate how this happens, consider, for simplicity, a non-rotating 
fluid (u0 = 0) in the Schwarzschild space-time. It may be easily shown that the acoustic metric takes the form 



ds ^ gtf 



1-1*2 



„ 1 
-dr - 2m- 



,2„,2 



1-1*2 



rdrdt 



1 2-c> 
gtt 1-1^2 



■dr 



2 ' rh 



(36) 



where gu = —(1 — 2/r), u — \vr\/ y^—gtt, and we have omitted the irrelevant conformal factor p/{hcs). Using the 
coordinate transformation 



dt^dt 



u 1 



9tt c% 



we remove the off-diagonal part from ([5^ and obtain 



ds' 



9tf 



1-1*2 



-df" 
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2 - C2li2 



1.2(1 -c2) 



.2\2 



dr^ + 



gtt L 1 - ""2 (c2 - u2)(l - u2)_ 
Next, we evaluate the metric near the acoustic horizon at r rh using the expansion in r — rh at first order 

d 



and making the substitution 



rh = 



2c.^ (cg 
or 



-gtt d_ 
2cs(l-c2) dr 



[r - Hi) 



(cs - u) 



(37) 



(38) 



(39) 



(40) 



where R denotes a new radial variable. Neglecting the first term in the square brackets in (|38p and setting t — ir, we 
obtain the Euclidean metric in the form 



ds| = n'^R^dr'^ + dR^ 



(41) 



where 



'9tt 



(42) 



Hence, the metric near r = rh is the product of the metric on and the Euclidean Rindler space-time 

ds\ = dR^ + R^dinrf. (43) 

With the periodic identification t = t + 27t/k, the metric (|43l) describes in plane polar coordinates. 

Furthermore, making the substitutions R = e'^^ / n and (p — y / r^ -\- t: , the Euclidean action takes the form of the 
2-f 1-dimensional free scalar field action at non-zero temperature 

5eM = / dT dx dy-{d,,^)\ (44) 

where we have set the upper and lower bounds of the integral over dy to -|-oo and — oo, respectively, assuming that rh 
is sufhciently large. Hence, the functional integral Z in (|35p is evaluated over the fields >p{x, y, r) that are periodic in 
T with period 27r/K. In this way, the functional Z is just the partition function for a grand-canonical ensemble of free 
bosons at the Hawking temperature Th = k/ {2ttkb)- However, the radiation spectrum will not be exactly thermal 
since we have to cut off the scales below the atomic scale (Unruh 1995). The choice of the cutoff and the deviation of 
the acoustic radiation spectrum from the thermal spectrum is closely related to the so-called transplanckian problem 
of Hawking radiation (Jacobson 1999a, 1992; Corley & Jacobson 1996). 

In the Newtonian approximation, (|42p reduces to the usual non-relativistic expression for the acoustic surface 
gravity represented by (pB]) . 
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IX. SALIENT FEATURES OF ACOUSTIC BLACK HOLES AND ITS CONNECTION TO 

ASTROPHYSICS 

In summary, analogue (acoustic) black holes (or systems) are fluid-dynamic analogue of general relativistic black holes. 
Analogue black holes possess analogue (acoustic) event horizons at local transonic points. Analogue black holes emit 
analogue Hawking radiation, the temperature of which is termed as analogue Hawking temperature, which may be 
computed using Newtonian description of fluid flow. Black hole analogues are important to study because it may be 
possible to create them experimentally in laboratories to study some properties of the black hole event horizon, and 
to study the experimental manifestation of Hawking radiation. 

According to the discussion presented in previous sections, it is now obvious that, to calculate the analogue 
surface gravity k and the analogue Hawking temperature Tj^fj for a classical analogue gravity system, one does 
need to know the exact location (the radial length scale) of the acoustic horizon r^, the dynamical and the acoustic 
velocity corresponding to the flowing fluid at the acoustic horizon, and its space derivatives, respectively. Hence an 
astrophysical fluid system, for which the above mentioned quantities can be calculated, can be shown to represent an 
classical analogue gravity model. 

For acoustic black holes, in general, the ergo-sphere and the acoustic horizon do not coincide. However, for some 
specific stationary geometry they do. This is the case, e.g. in the following two examples: 

1. Stationary spherically symmetric configuration where fluid is radially falling into a pointlike drain at the origin. 
Since u = u±_ everywhere, there will be no distinction between the ergo-sphere and the acoustic horizon. An 
astrophysical example of such a situation is the stationary spherically symmetric Bondi-type accretion (Bondi 
1952) onto a Schwarzschild black hole, or onto other non rotating compact astrophysical objects in general, see 
section 10.2 for further details on spherically symmetric astrophysical accretion. 

2. Two-dimensional axisymmctric configuration, where the fluid is radially moving towards a drain placed at the 
origin. Since only the radial component of the velocity is non-zero, u = u± everywhere. Hence, for this system, 
the acoustic horizon will coincide with the ergo region. An astrophysical example is an axially symmetric 
accretion with zero angular momentum onto a Schwarzschild black hole or onto a non-rotating neutron star, see 
section 10.3 for further details of axisymmctric accretion. 

In subsequent sections, we thus concentrate on transonic black hole accretion in astrophysics. We will first review 
various kind of astrophysical accretion, emphasizing mostly on the black hole accretion processes. We will then show 
that sonic points may form in such accretion and the sonic surface is essentially an acoustic horizon. We will provide 
the formalism using which one can calculate the exact location of the acoustic horizon (sonic points) r/j, the dynamical 
accretion velocity u and the acoustic velocity Cg at r^, and the space gradient of those velocities (du/dr) and {dcg/dr) 
at rh, respectively. Using those quantities, we will then calculate k and Tah for an accreting black hole system. 
Such calculation will ensure that accretion processes in astrophysics can be regarded as a natural example of classical 
analogue gravity model. 

X. TRANSONIC BLACK HOLE ACCRETION IN ASTROPHYSICS 

A. A General Overview 

Gravitational capture of surrounding fluid by massive astrophysical objects is known as accretion. There remains a 
major difference between black hole accretion and accretion onto other cosmic objects including neutron stars and 
white dwarfs. For celestial bodies other than black holes, infall of matter terminates either by a direct collision with 
the hard surface of the accretor or with the outer boundary of the magneto-sphere, resulting the luminosity (through 
energy release) from the surface. Whereas for black hole accretion, matter ultimately dives through the event horizon 
from where radiation is prohibited to escape according to the rule of classical general relativity, and the emergence 
of luminosity occurs on the way towards the black hole event horizon. The efficiency of accretion process may be 
thought as a measure of the fractional conversion of gravitational binding energy of matter to the emergent radiation, 
and is considerably high for black hole accretion compared to accretion onto any other astrophysical objects. Hence 
accretion onto classical astrophysical black holes has been recognized as a fundamental phenomena of increasing 
importance in relativistic and high energy astrophysics. The extraction of gravitational energy from the black hole 
accretion is believed to power the energy generation mechanism of X-ray binaries and of the most luminous objects of 
the Universe, the Quasars and active galactic nuclei (Frank, King & Raine 1992). The black hole accretion is, thus, 
the most appealing way through which the all pervading power of gravity is explicitly manifested. 
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As it is absolutely impossible to provide a detail discussion of a topic as vast and diverse as accretion onto various 
astrophysical objects in such a small span, this section will mention only a few topic and will concentrate on fewer 
still, related mostly to accretion onto black hole. For details of various aspects of accretion processes onto compact 
objects, recent reviews like Pringle 1981; Chakrabarti 1996a; Wiita 1998; Lin & Papaloizou 1996; Blandford 1999; 
Rees 1997; Bisnovayati-Kogan 1998; Abramowicz et al 1998; and the monographs by Frank, King & Raine 1992, and 
Kato, Fukue & Mineshige 1998, will be of great help. 

Accretion processes onto black holes may be broadly classified into two different categories. When accreting 
material does not have any intrinsic angular momentum, flow is spherically symmetric and any parameters governing 
the accretion will be a functio n of radial distance only. On the other hand, for matter accreting with considerable 
intrinsic angular momentum, [226j | flow geometry is not that trivial. In this situation, before the infalling matter 
plunges through the event horizon, accreting fluid will be thrown into circular orbits around the hole, moving inward 
usually when viscous stress in the fluid helps to transport away the excess amount of angular momentum. This outward 
viscous transport of angular momentum of the accreting matter leads to the formation of accretion disc around the 
hole. The structure and radiation spectrum of these discs depends on various physical parameters governing the flow 
and on specific boundary conditions. 

If the instantaneous dynamical velocity and local acoustic velocity of the accreting fluid, moving along a space 
curve parameterized by r, are u(r) and Cs{r), respectively, then the local Mach number M{r) of the fluid can be 
defined as M{r) — u{r)/cs{r). The flow will be locally subsonic or supersonic according to M{r) < 1 or > 1, i.e., 
according to uir) < Cs(r) or u(r) > Cs(r). The flow is transonic if at any moment it crosses M — 1. This happens when 
a subsonic to supersonic or supersonic to subsonic transition takes place either continuously or discontinuously. The 
point (s) where such crossing takes place continuously is (are) called sonic point (s), and where such crossing takes place 
discontinuously are called shocks or discontinuities. At a distance far away from the black hole, accreting material 
almost always remains subsonic (except for the supersonic stellar wind fed accretion) since it possesses negligible 
dynamical flow velocity. On the other hand, the flow velocity will approach the velocity of light (c) while crossing 
the event horizon, while the maximum possible value of sound speed (even for the steepest possible equation of state) 
would be c/-\/3, resulting M > 1 close to the event horizon. In order to satisfy such inner boundary condition imposed 
by the event horizon, accretion onto black holes exhibit transonic properties in general. 



M<1 KcO^®*'*^ '^ST/. 



1>1 



M<1 



FIG. 1: Spherically symmetric transonic black hole accretion with acoustic horizon. 
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B. Mono-transonic Spherical Accretion 

Investigation of accretion processes onto celestial objects was initiated by Hoyle & Lyttleton (1939) by computing 
the rate at which prcssurc-lcss matter would be captured by a moving star. Subsequently, theory of stationary, 
spherically symmetric and transonic hydrodynamic accretion of adiabatic fluid onto a gravitating astrophysical object 
at rest was formulated in a seminal paper by Bondi (1952) using purely Newtonian potential and by including 
the pressure effect of the accreting material. Later on, Michel (1972) discussed fully general relativistic polytropic 
accretion on to a Schwarzschild black hole by formulating the governing equations for steady spherical flow of perfect 
fluid in Schwarzschild metric. FoUowing Michel's relativistic generalization of Bondi's treatment, Bcgelman (1978) 
and Moncrief (1980) discussed some aspects of the sonic points of the flow for such an accretion. Spherical accretion 
and wind in general relativity have also been considered using equations of state other than the polytropic one 
and by incorporating various radiative processes (Shapiro 1973, 1973a; Blumenthal & Mathews 1976; Brinkmann 
1980). Malec (1999) provided the solution for general relativistic spherical accretion with and without back reaction, 
and showed that relativistic effects enhance mass accretion when back reaction is neglected. The exact values of 
dynamical and thermodynamic accretion variables on the sonic siirfacc, and at extreme close vicinity of the black hole 
event horizons, have recently been calculated using complete general relativistic (Das 2002) as well as pseudo general 
relativistic (Das & Sarkar 2001) treatments. 

Figure 1 pictorially illustrates the generation of the acoustic horizon for spherical transonic accretion. Let us 
assume that an isolated black hole at rest accretes matter. The black hole (denoted by B in the figure) is assumed to 
be of Schwarzschild type, and is embedded by an gravitational event horizon of radius 2GMbh / c? . Infalling matter 
is assumed not to possess any intrinsic angular momentum, and hence, falls freely on to the black hole radially. 
Such an accreting system possesses spherical symmetry. Far away from the black hole the dynamical fluid velocity 
is negligible and hence the matter is subsonic, which is demonstrated in the figure by M < 1. In course of its 
motion toward the event horizon, accreting material acquires sufficiently large dynamical velocity due to the black 
hole's strong gravitational attraction. Consequently, at a certain radial distance, the Mach number becomes unity. 
The particular value of r, for which M=l, is referred as the transonic point or the sonic point, and is denoted by 
Thi as mentioned in the above section. For r < rh, matter becomes supersonic and any acoustic signal created in 
that region is bound to be dragged toward the black hole, and can not escape to the region r > r^- In other words, 
any co-moving observer from r<r/i can not communicate with another observer at r > by sending any signal 
traveling with velocity Wsignai<Cs. Hence the spherical surface through rn is actually an acoustic horizon for stationary 
configuration, which is generated when accreting fluid makes a transition from subsonic (M < 1) to the supersonic (M 
> 1) state. In subsequent sections, we will demonstrate how one can determine the location of rh and how the surface 
gravity and the analogue Hawking temperature corresponding to such r/, can be computed. Note, however, that for 
spherically symmetric accretion, only one acoustic horizon may form for a given set of initial boundary configuration 
characterizing the stationary configuration. For matter accreting with non-zero intrinsic angular momentum, multiple 
acoustic horizons can be obtained. Details of such configurations will be discussed in subsequent sections. 

It is perhaps relevant to mention that spherical black hole accretion can allow standing shock formation. Pertur- 
bations of various kinds may produce discontinuities in an astrophysical fluid flow. By discontinuity at a surface in a 
fluid flow we understand any discontinuous change of a dynamical or a thermodynamic quantity across the surface. 
The corresponding surface is called a surface of discontinuity. Certain boundary conditions must be satisfied across 
such surfaces and according to these conditions, surfaces of discontinuities are classified into various categories. The 
most important such discontinuities are shock waves or shocks. 

While the possibility of the formation of a standing spherical shock around compact objects was first conceived 
long ago (Bisnovatyi-Kogan, Zel'Dovich, & Sunyaev 1971), most of the works on shock formation in spherical accretion 
share more or less the same philosophy that one should incorporate shock formation to increase the efficiency of directed 
radial infall in order to explain the high luminosity of AGNs and QSOs and to model their broad band spectrum 
(Jones & Ellison 1991). Considerable work has been done in this direction where several authors have investigated the 
formation and dynamics of standing shock in spherical accretion (Mcszaros & Ostriker 1983; Protheros & Kazanas 
1983; Chang & Osttriker 1985; Kazanas & Ellision 1986; Babul, Ostriker & Meszaros 1989; Park 1990, 1990a). 

Study of spherically symmetric black hole accretion leads to the discovery of related interesting problems like 
entropic-acoustic or various other instabilities in spherical accretion (Foglizzo & Tagger 2000; Blondin & Ellison 2001; 
Lai & Goldreich 2000; Foglizzo 2001; Kovalenko & Eremin 1998), the realizability and the stability properties of 
Bondi solutions (Ray & Bhattacharjee 2002), production of high energy cosmic rays from AGNs (Protheroe & Szabo 
1992), study of the hadronic model of AGNs (Blondin & Konigl 1987; Contopoulos & Kazanas 1995), high energetic 
emission from relativistic particles in our galactic centre (Markoff, Melia & Sarcevic 1999), explanation of high lithium 
abundances in the late-type, low-mass companions of the soft X-ray transient, (Guessomn & Kazanas 1999), study of 
accretion powered spherical winds emanating from galactic and extra galactic black hole environments (Das 2001). 
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C. Breaking the Spherical Symmetry: Accretion Disc 

1. A General Overview 

In sixties, possible disc-like structures around one of the binary components were found (Kraft, 1963) and some 
tentative suggestions that matter should accrete in the form of discs were put forward (Pcndergest & Burbidge 1968; 
Lynden-Bell 1969). Meanwhile, it was understood that for spherically symmetric accretion discussed above, the 
(radial) infall velocity is very high, hence emission from such a rapidly falling matter was not found to be strong 
enough to explain the high luminosity of Quasars and AGNs. Introducing the idea of magnetic dissipation, efforts 
were made to improve the luminosity (Shvartsman 1971, 1971a; Shapiro 1973, 1973a). 

Theoretically, accretion discs around black holes were first envisaged to occur within a binary stellar system where 
one of the components is compact object (i.e., white dwarfs, neutron stars or a black hole) and the secondary would 
feed matter onto the primary either through an wind or through Roche lobe overflow. In either case, the accreted 
matter would clearly possesses substantial intrinsic angular momentum with respect to the compact object (a black 
hole, for our discussion). A flow with that much angular momentum will have much smaller infall velocity and much 
higher density compared to the spherical accretion. The infall time being higher, viscosity within the fluid, presumably 
produced by turbulence or magnetic field, would have time to dissipate angular momentum (except in regions close 
to the black holes, since large radial velocity close the event horizon leads to the typical value of dynamical time 
scale much smaller compared to the viscous time scale) and energy. As matter loses angular momentum, it sinks 
deeper into the gravitational potential well and radiate more efficiently. The flow encircles the compact accretor and 
forms a quasi-stationary disc like structure around the compact object and preferably in the orbital plane of it. Clear 
evidences for such accretion discs around white dwarfs in binaries was provided by analysis of Cataclysmic variable 
(Robinson 1976). 

Accretion forming a Keplarian disc [227| around a Schwarzschild black hole produces efficiency r/ (the fraction 
of gravitational energy released) of the order of 77 '--^ 0.057 and accretion onto a maximally rotating Kerr black hole 
is even more efficient, yielding rj ~ 0.42. However, the actual efficiencies depends on quantities such as viscosity 
parameters and the cooling process inside the disc (see Wiita 1998 and references therein). This energy is released in 
the entire electromagnetic spectrum and the success of a disc model depends on its ability to describe the way this 
energy is distributed in various frequency band. 

In case of binary systems, where one of the components is a compact object like white dwarfs, neutron star 
or a black hole, the companion is stripped off its matter due to the tidal effects. The stripped off matter, with 
angular momentum equal to that of the companion, gradually falls towards the central compact object as the angular 
momentum is removed by viscosity. As the flow possesses a considerable angular momentum to begin with, it is 
reasonable to assume that the disc will form and the viscosity would transport angular momentum from inner part 
of the disc radially to the outer part which allows matter to further fall onto the compact body. This situation could 
be described properly by standard thin accretion disc, which may be Keplarian in nature. On the other hand, in the 
case of active galaxies and quasars, the situation could be somewhat different. The supermassive {Mbh ^ 10^ Mq) 
central black hole is immersed in the intergalactic matter. In absence of any binary companion, matter is supplied to 
the central black hole very intermittently, and the angular momentum of the accreting matter at the outer edge of the 
disc may be sub-Keplarian. This low angular momentum flow departs the disc from Keplarian in nature and a 'thick 
disc' is more appropriate to describe the behaviour instead of standard thin, Keplarian Shakura Sunyaev (Shakura & 
Sunyaev 1973) disc. 

2. Thm Disc Model 

In standard thin disc model (Shakura & Sunyaev 1973; Novikov & Thorne 1973), originally conceived to describe 
Roche lobe accretion in a binary system, the local height H{r) of the disc is assumed to be small enough compared to 
the local radius of the disc r, i.e., the 'thinness' condition is dictated by the fact that H{r) « r. Pressure is neglected 
so that the radial force balance equations dictates the specific angular momentum distribution to become Keplarian 
and the radial velocity is negligible compared to the azimuthal velocity (u^ << Vrj,). Unlike the spherical accretion, 
temperature distribution is far below than virial. Under the above mentioned set of assumptions, radial equations 
of steady state disc structure could be decoupled from the vertical ones and could be solved independently. The 
complete solutions describing the steady state disc structure can be obtained by solving four relativistic conservation 
equations, namely; the conservation of rest mass, specific angular momentum, specific energy and vertical momentum 
balance condition. In addition, a viscosity law may be specified which may transport angular momentum outwards 
allowing matter to fall in. On the top of it, in standard thin disc model, the shear is approximated as proportional to 
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the pressure of the disc with proportionahty constant a, a being the viscosity parameter having numerical value less 
than unity. 

High uncertainty remains in investigating the exact nature of the viscosity inside a thin accretion disc (see Wiita 
1998 and references therein). One of the major problems is to explain the origin of sufficiently large viscosity that 
seems to be present inside accretion discs in the binary system. Unfortunately, under nearly all astrophysically 
relevant circumstances, all of the well understood microscopic transverse momentum transport mechanism such as 
ionic, molecular and radiative viscosity are extremely small. Observations with direct relevance to the nature and 
strength of the macroscopic viscosity mechanism are very difficult to make; the only fairly direct observational evidence 
for the strength of disc viscosity comes from the dwarf novae system. For a black hole as compact accretor, such 
observational evidences is far from reality till date. Therefore advances in understanding the disc viscosity is largely 
based on theoretical analysis and numerical techniques. Usually accepted view is that the viscosity may be due to 
magnetic transport of angular momentum or due to small scale turbulent dissipation. Over the past several years an 
explanation of viscosity in terms of Velikhov-Chandrasekhar-Balbus-Hawley instability (linear magnetic instability) 
has been investigated; see, e.g., Balbus & Hawle 1998 for further details. 

3. Thick Disc Model 

The assumptions implying accretion discs are always thin can break down in the innermost region. Careful consider- 
ation of the effects of general relativity show that the flow must go supersonically through a cusp. For considerably 
high accretion rate, radiation emitted by the in-falling matter exerts a significant pressure on the gas. The radiation 
pressure inflates the disc, and make it geometrically thick {H{r) ~ r, at least for the inner 10 — lOOrg), which is often 
otherwise known as 'accretion torus'. This considerable amount of radiation pressure must be incorporated to find the 
dynamical structure of the disc and in determining the thermodynamical quantities inside the disc. Incorporation of 
the radiation pressure term in Euler equation dictates the angular momentum deviation from that of the Keplarian. 
The angular momentum distribution becomes super (sub) Keplarian if the pressure gradient is positive (negative). 

Introducing a post-Newtonian (these 'pseudo' potentials are widely used to mimic the space time around the 
Schwarzschild or the Kerr metric very nicely, see section 14 for details) $ = — ^^2r') '^^ ^^^^ °^ usual $ Newtonian = 
— QMsjL (where Vg is the 'gravitational' radius), Paczynski and Wiita (1980) provided the first thick disc model which 
joins with the standard thin disc at large radius without any discontinuity. They pointed out several important features 
of these configuration. It has been shown that the structure of thick disc in inner region is nearly independent of 
the viscosity and efficiency of accretion drops dramatically. More sophisticated model of radiation supported thick 
disc including self-gravity of the disc with full general relativistic treatment was introduced later (Wiita 1982; Lanza 
1992). 

4.. Further Developments 

Despite having a couple of interesting features, standard thick accretion disc model suffers from some limitations for 
which its study fell from favour in the late '80s. Firstly, the strong anisotropic nature of the emission properties of 
the disc has been a major disadvantage. Secondly, a non-accreting thick disc is found to be dynamically and globally 
unstable to non-axisymmetric perturbations. However, an ideal 'classical thick disc', if modified to incorporate high 
accretion rates involving both low angular momentum and considerable radial infall velocity self-consistently, may 
remain viable. Also, it had been realized that neither the Bondi (1952) fiow nor the standard thin disc model 
could individually fit the bill completely. Accretion disc theorists were convinced about the necessity of having an 
intermediate model which could bridge the gap between purely spherical flow (Bondi type) and purely rotating flow 
(standard thin disc). Such modification could be accomplished by incorporating a self-consistent 'advection' term 
which c:oul(l take care of finite radial velocity of accreting material (for the black hole candidates wliic;li may gradually 
approaches the velocity of light to satisfy the inner boundary condition on event horizon) along with its rotational 
velocity and generalized heating and cooling terms (Hoshi fc Shibazaki 1977; Liang & Thompson 1980; Ichimaru 1977; 
Paczynski & Bisnobhatyi-Kogan 1981; Abramowicz & Zurck 1981; Muchotrzeb & Paczynski 1982; Muchotrzeb 1983; 
Fukue 1987; Abramowicz et al. 1988; Narayan & Yi 1994; Chakrabarti 1989, 1996). 

D. Multi-transonic Accretion Disc 

For certain values of the intrinsic angular momentum density of accreting material, the number of sonic point, 
unlike spherical accretion, may exceed one, and accretion is called 'multi-transonic'. Study of such multi-transonicity 
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was initiated by Abramowicz & Zurck (1981). Subsequently, multi-transonic accretion disc has been studied in a 
number of works (Fukue 1987; Cliakrabarti 1990, 1996; Kafatos & Yang 1994; Yang & Kafatos 1995; Pariev 1996; 
Peitz & Appl 1997; Lasota & Abramowicz 1997; Lu, Yu, Yuan & Young 1997; Das 2004; Barai, Das & Wiita 2004; 
Abraham, Bilic & Das 2006; Das, Bilic & Dasgupta 2006). All the above works, except Barai, Das & Wiita 2004, 
usually deal with low angular momentum sub-Keplerian inviscid flow around a Schwarzschild black hole or a prograde 
flow around a Kerr black hole. Barai, Das & Wiita 2004 studied the retrograde flows as well and showed that a higher 
angular momentum (as high as Keplerian) retrograde flow can also produce multi-transonicity. Sub-Keplerian weakly 
rotating flows are exhibited in various physical situations, such as detached binary systems fed by accretion from 
OB stellar winds (lUarionov & Sunyacv 1975; Liang & Nolan 1984), scmi-dctachcd low-mass non-magnetic binaries 
(Bisikalo et al. 1998), and super-massive black holes fed by accretion from slowly rotating central stellar clusters 
(lUarionov 1988; Ho 1999 and references therein). Even for a standard Keplerian accretion disc, turbulence may 
produce such low angular momentum flow (see, e.g., Igumenshchev & Abramowicz 1999, and references therein). 



E. Non-axisymmetric Accretion Disc 

All the above mentioned works deals with 'axisymmetric' accretion, for which the orbital angular momentum of 
the entire disc plane remains aligned with the spin angular momentum of the compact object of our consideration. 
In a strongly coupled binary system (with a compact object as one of the components), accretion may experience a 
non-axisymmetric potential because the secondary donor star may exert non-axisymmetric tidal force on the accretion 
disc around the compact primary. In general, non-axisymmetric tilted disc may form if the accretion takes place out 
of the symmetry plane of the spinning compact object. Matter in such misaligned disc will experience a torque due 
to the general relativistic Lcnsc-Thirring effect (Lense & Thirring 1918), leading to the precession of the inner disc 
plane. The differential precession with radius may cause stress and dissipative effects in the disc. If the torque 
remains strong enough compared to the internal viscous force, the inner region of the initially tilted disc may be 
forced to realigned itself with the spin angular momentum (symmetry plane) of the central accretor. This phenomena 
of partial re-alignment (out to a certain radial distance known as the 'transition radius' or the 'alignment radius') of 
the initially non-axisymmetric disc is known as the 'Bardeen-Petterson effect' (Bardeen & Petterson 1975). Such a 
transition radius can be obtained by balancing the precession and the inward drift or the viscous time scale. 

Astrophysical accretion disc subjected to the Bardeen-Petterson effect becomes 'twisted' or 'warped'. A large 
scale warp (twist) in the disc may modify the emergent spectrum and can influence the direction of the Quasar and 
micro-quasar jets emanating out from the inner region of the accretion disc (see, e.g., Maccarone 2002; Lu & Zhou 
2005, and references therein). 

Such a twisted disc may be thought as an ensemble of annuli of increasing radii, for which the variation of the 
direction of the orbital angular momentum occurs smoothly while crossing the alignment radius. System of equations 
describing such twisted disc have been formulated by several authors (see, e.g., Peterson 1977; Kumar 1988; Demianski 
& Ivanov 1997; and references therein), and the time scale required for a Kerr black hole to align its spin angular 
momentum with that of the initially misaligned accretion disc, has also been estimated (Scheuer & Feiler 1996). 
Numerical simulation using throe dimensional Newtonian Smooth Particle Hydrodynamics (SPH) code (Nelson & 
Papaloizou 2000) as well as using fully general relativistic framework (Fragile & Anninos 2005) reveal the geometric 
structure of such discs. 

We would, however, not like to explore the non-axisymmetric accretion further in this review. One of the main 
reasons for which is, as long as the acoustic horizon forms at a radial length scale smaller than that of the alignment 
radius (typically 100 Vg - 1000 rg, according to the original estimation of Bardeen & Petterson 1975), one need not 
implement the non-axisymmetric geometry to study the analogue effects. 



F. Angular Momentum Supported Shock in Multi-transonic Accretion Disc 

In an adiabatic flow of the Newtonian fluid, the shocks obey the following conditions (Landau & Lifshitz 1959) 



[[H] = o, [[p + pu^]]=o, 



0, 



where [[/]] denotes the discontinuity of / across the surface of discontinuity, i.e. 

[[/]] = /2 - /l, 



(45) 
(46) 



with /2 and /i being the boundary values of the quantity / on the two sides of the surface. Such shock waves arc quite 
often generated in various kinds of supersonic astrophysical flows having intrinsic angular momentum, resulting in a 
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flow which becomes subsonic. This is because the repulsive centrifugal potential barrier experienced by such flows is 
sufflciently strong to brake the infalling motion and a stationary solution could be introduced only through a shock. 
Rotating, transonic astrophysical fluid flows are thus believed to be 'prone' to the shock formation phenomena. 

One also expects that a shock formation in black-hole accretion discs might be a general phenomenon because 
shock waves in rotating astrophysical flows potentially provide an important and efficient mechanism for conversion 
of a significant amount of the gravitational energy into radiation by randomizing the directed infall motion of the 
accreting fluid. Hence, the shocks play an important role in governing the overall dynamical and radiative processes 
taking place in astrophysical fluids and plasma accreting onto black holes. The study of steady, standing, stationary 
shock waves produced in black hole accretion has acquired an important status, and a number of works studied the 
shock formation in black hole accretion discs (Fukue 1983; Hawley, Wilson & Smarr 1984; Ferrari et al. 1985; Sawada, 
Matsuda & Hachisu 1986; Spruit 1987; Chakrabarti 1989; Abramowicz & Chakrabarti 1990; Yang & Kafatos 1995; 
Chakrabarti 1996a; Lu, Yu, Yuan & Young 1997; Caditz & Tsuruta 1998; Toth, Keppens & Botchev 1998; Das 2002; 
Takahashi, Rillet, Fukumura & Tsuruta 2002; Das, Pendharkar & Mitra 2003; Das 2004; Chakrabarti & Das 2004; 
Fukumura & Tsuruta 2004; Abraham, Bilic & Das 2006; Das, Bilic & Dasgupta 2006) For more details and for a more 
exhaustive list of references see, e.g., Chakrabarti 1996c and Das 2002. 

Generally, the issue of the formation of steady, standing shock waves in black-hole accretion discs is addressed 
in two different ways. First, one can study the formation of Rankine-Hugoniot shock waves in a polytropic flow. 
Radiative cooling in this type of shock is quite inefficient. No energy is dissipated at the shock and the total specific 
energy of the accreting material is a shock-conserved quantity. Entropy is generated at the shock and the post-shock 
flow possesses a higher entropy accretion rate than its pre-shock counterpart. The flow changes its temperature 
permanently at the shock. Higher post-shock temperature puffs up the post-shock flow and a quasi-spherical, quasi- 
toroidal centrifugal pressure supported region is formed in the inner region of the accretion disc (see Das 2002, and 
references therein for further detail) which locally mimics a thick accretion flow. 

Another class of the shock studies concentrates on the shock formation in isothermal black-hole accretion discs. 
The characteristic features of such shocks are quite different from the non-dissipative shocks discussed above. In 
isothermal shocks, the accretion flow dissipates a part of its energy and entropy at the shock surface to keep the post- 
shock temperature equal to its pre-shock value. This maintains the vertical thickness of the flow exactly the same just 
before and just after the shock is formed. Simultaneous jumps in energy and entropy join the pre-shock supersonic 
flow to its post-shock subsonic counterpart. For detailed discussion and references see, e.g., Das, Pendharkar & Mitra 
2003, and Fukumura & Tsuruta 2004. 

In section 13.5, wc will construct and solve the equations governing the general relativistic Rankine-Hugoniot 
shock. The shocked accretion flow in general relativity and in post-Newtonian pseudo-Schwarzschild potentials will 
be discussed in the section 13.5 - 13.8 and 16.2 respectively. 

XI. MOTIVATION TO STUDY THE ANALOGUE BEHAVIOUR OF TRANSONIC BLACK HOLE 

ACCRETION 

Since the publication of the seminal paper by Bondi in 1952 (Bondi 1952), the transonic behaviour of accreting fluid 
onto compact astrophysical objects has been extensively studied in the astrophysics community, and the pioneering 
work by Unruh in 1981 (Unruh 1981), initiated a substantial number of works in the theory of analogue Hawking 
radiation with diverse fields of application stated in section 4-5. It is surprising that no attempt was made to bridge 
these two categories of research, astrophysical black hole accretion and the theory of analogue Hawking radiation, by 
providing a self-consistent study of analogue Hawking radiation for real astrophysical fluid flows, i.e., by establishing 
the fact that accreting black holes can be considered as a natural example of analogue system. Since both the theory 
of transonic astrophysical accretion and the theory of analogue Hawking radiation stem from almost exactly the 
same physics, the propagation of a transonic fluid with acoustic disturbances embedded into it, it is important to 
study analogue Hawking radiation for transonic accretion onto astrophysical black holes and to compute Tah for such 
accretion. 

In the following sections, we will describe the details of the transonic accretion and will show how the accreting 
black hole system can be considered as a classical analogue system. We will first discuss general relativistic accretion of 
spherically symmetric (mono-transonic Bondi (1952) type accretion) and axisymmetric (multi-transonic disc accretion) 
flow. We will then introduce a number of post-Newtonian pseudo-Schwarzschild black hole potential, and will discuss 
black hole accretion under the influence of such modified potentials. 
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XII. GENERAL RELATIVISTIC SPHERICAL ACCRETION AS AN ANALOGUE GRAVITY MODEL 

In this section, we will demonstrate how one can construct and solve the equations governing the general relativistic, 
spherically symmetric, steady state accretion flow onto a Schwarzschild black hole. This section is largely based on 
Das 2004a. 

Accretion flow described in this section is 9 and (/> symmetric and possesses only radial inflow velocity. In this 
section, we use the gravitational radius Vg as rg — 2GMbh /c^ ■ The radial distances and velocities are scaled in 
units of rg and c respectively and all other derived quantities are scaled accordingly; G = c = Mbh = 1 is used. 
Accretion is governed by the radial part of the general relativistic time independent Euler and continuity equations 
in Schwarzschild metric. We will consider the stationary solutions. We assume the dynamical in-fall time scale to be 
short compared with any dissipation time scale during the accretion process. 



A. The Governing Equations 



To describe the fluid, we use a polytropic equation of state (this is common in the theory of relativistic black hole 
accretion) of the form 



(47) 



where the polytropic index 7 (equal to the ratio of the two specific heats Cp and c^) of the accreting material is 
assumed to be constant throughout the fluid. A more realistic model of the flow would perhaps require a variable 
polytropic index having a functional dependence on the radial distance, i.e. 7=7(r). However, we have performed 
the calculations for a sufRciently large range of 7 and we believe that all astrophysically relevant polytropic indices 
are covered in our analysis. 

The constant K in ((47|) may be related to the specific entropy of the fluid, provided there is no entropy generation 
during the flow. If in addition to (jiT]) the Clapeyron equation for an ideal gas holds 



P 



Kb 
firUp' 



(48) 



where T is the locally measured temperature, fi the mean molecular weight, mn^rnp the mass of the hydrogen atom, 
then the specific entropy, i.e. the entropy per particle, is given by (Landau & Lifshitz 1959): 



1 



7-1 



\ogK 



7 



7-1 



+ constant. 



(49) 



where the constant depends on the chemical composition of the accreting material. Equation (j49p confirms that K 
in (|47p is a measure of the specific entropy of the accreting matter. 

The specific enthalpy of the accreting matter can now be defined as 



h 



(p + e) 



(50) 



where the energy density e includes the rest-mass density and the internal energy and may be written as 

P 



The adiabatic speed of sound is defined by 



From (ISTI) we obtain 



7-1 



dp 



dp 



constant entropy 



7 - 1 - 



7-1 



(51) 



(52) 



(53) 



Combination of ([5^ and (ITf|) gives 



(54) 
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Using the above relations, one obtains the expression for the specific enthalpy 

7-1 



h = 



7 - 1 - c2 



(55) 



The rest-mass density p, the pressure p, the temperature T of the flow and the energy density e may be expressed in 
terms of the speed of sound as 



p = K- 



7 - 1 - c2 



(56) 



p = K ^'-^ 



7 - 1 - c2 



(57) 



T : 



KB 

/imp 



7 - 1 - c2 



(58) 



1 /7-l\^-i 
\ - — ^ 

7 



7-1-c^ 



1 

1 + - 



7 V7- 1 



(59) 



The conserved specific fiow energy £ (the relativistic analogue of Bernoulli's constant) along each stream line 
reads £ = hut, (Anderson 1989) where h and are the specific enthalpy and the four velocity, which can be re-cast 
in terms of the radial three velocity u and the polytropic sound speed Cg to obtain: 



E = 



7-1 



7-(l + c2) 



1 - 1/? 



1 



(60) 



One concentrates on positive Bernoulli constant solutions. The mass accretion rate M may be obtained by integrating 
the continuity equation: 



r- 1 

■(1-U2) 



(61) 



where p is the proper mass density. 

We define the 'entropy accretion rate' S as a quasi-constant multiple of the mass accretion rate in the following 
way: 



1 



S = if 1 - 7 M = inpur^ 



I r- 1 



cl (7-1) 



7-(l + c2) 



(62) 



Note that, in the absence of creation or annihilation of matter, the mass accretion rate is a universal constant of 
motion, whereas the entropy accretion rate is not. As the expression for S contains the quantity K = p/p'^, which 
measures the specific entropy of the fiow, the entropy rate S remains constant throughout the flow only if the entropy 
per particle remains locally unchanged. This latter condition may be violated if the accretion is accompanied by a 
shock. Thus, S is a constant of motion for shock-free polytropic accretion and becomes discontinuous (increases) at 
the shock location, if a shock forms in the accretion. One can solve the two conservation equations for £ and S to 
obtain the complete accretion profile. 



B. Transonicity 

Simultaneous solution of (|60ll62p provides the dynamical three velocity gradient at any radial distance r: 

(l-u^) [c2(4r-3)-l] _Afir,u,Cs) 



du 
dr 



2r (r-l)(M2-c2) 



(63) 
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A real physical transonic flow must be smooth everywhere, except possibly at a shock. Hence, if the denominator 
T) (r, u, Cg) of (j63p vanishes at a point, the numerator Af (r, u, c^) must also vanish at that point to ensure the physical 
continuity of the flow. Borrowing the terminology from the dynamical systems theory (see, e.g., Jordan & Smith 2005), 
one therefore arrives at the critical point conditions by making V (r, u, Cg) and Af (r, u, Cs) of (|63p simultaneously equal 
to zero. We thus obtain the critical point conditions as: 



l(r= 



-■s\{r 



1 



(64) 



Tc being the location of the critical point or the so called 'fixed point' of the differential equation ([55)1 . 

From (|64p . one easily obtains that Mc, the Mach number at the critical point, is exactly equal to unity. This 
ensures that the critical points are actually the sonic points, and thus, rc is actually the location of the acoustic event 
horizon. In this section, hereafter, we will thus use Vh in place of rc- Note, however, that the equivalence of the 
critical point with the sonic point (and thus with the acoustic horizon) is not a generic feature. Such an equivalence 
strongly depends on the flow geometry and the equation of state used. For spherically symmetric accretion (using 
any equation of state) , or polytropic disc accretion where the expression for the disc height is taken to be constant 
(Abraham, Bilic & Das 2006), or isothermal disc accretion with variable disc height, such an equivalence holds good. 
For all other kind of disc accretion, critical points and the sonic points are not equivalent, and the acoustic horizon 
forms at the sonic points and not at the critical point. We will get back to this issue in greater detail in section 13.3. 



Substitution of u|(r=rh) ^^'^ Cs |(r=rh) i^i^o ()60p for r = rh provides: 

rlTi + rhT2 + Fg = 



(65) 



where 



£2 (37-2)' 



Fi = 
-27(7 



2£2 (2 „ 3^) + 9 (^ _ 1) 



1 



32 (£2 - 1) (7 - 1) 



4(7-1) (£2-1) 

2 



27 



64 (f 2 _ 1) 



(66) 



Solution of (j65|) provides the location of the acoustic horizon in terms of only two accretion parameters {£, 7}, which 
is the two parameter input set to study the flow. 

We now set the appropriate limits on {£ , 7}to model the realistic situations encountered in astrophysics. As £ is 
scaled in terms of the rest mass energy and includes the rest mass energy, £ < 1 corresponds to the negative energy 
accretion state where radiative extraction of rest mass energy from the fluid is required. For such extraction to be 
made possible, the accreting fluid has to possess viscosity or other dissipative mechanisms, which may violate the 
Lorentzian invariance. On the other hand, although almost any £ > 1 is mathematically allowed, large values of £ 
represents flows starting from infinity with extremely high thermal energy (see section 13.4 for further detail), and 
£ > 2 accretion represents enormously hot flow configurations at very large distance from the black hole, which are 
not properly conceivable in realistic astrophysical situations. Hence one sets 1 ^ £ ^ 2. Now, 7 = 1 corresponds 
to isothermal accretion where accreting fluid remains optically thin. This is the physical lower limit for 7, and 
7 < 1 is not realistic in accretion astrophysics. On the other hand, 7 > 2 is possible only for superdense matter 
with substantially large magnetic field (which requires the accreting material to be governed by general relativistic 
magneto-hydrodynamic equations, dealing with which is beyond the scope of this article) and direction dependent 
anisotropic pressure. One thus sets 1 ;S 7 ;S 2 as well, so {£,7}has the boundaries 1 ^ {£:l} ~ 2. However, one 
should note that the most preferred values of 7 for realistic black hole accretion ranges from 4/3 to 5/3 (Frank, King 
& Raine 1992). 

For any specific value of {£, 7}, (j65p can be solved completely analytically by employing the Cardano-Tartaglia-del 
Ferro technique. One defines: 



Si = 



3r2 - Tl 



S2 = 



gFiFa - 27F3 - 2r3 



f2i 
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Si 



S?, ^2 = {/S2- 



E^, e = cos- 



si 



(67) 



so that the three roots for rh come out to be: 



Fi 
3 
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However, note that not all V/i{i = 1, 2, 3} would be real for all {£, 7}. It is easy to show that if > 0, only one root 
is real; if \I> = 0, all roots are real and at least two of them are identical; and if < 0, all roots are real and distinct. 
Selection of the real physical {vh has to be greater than unity) roots requires a close look at the solution for rh for the 
astrophysically relevant range of {£,7}. One finds that for the preferred range of {£,7}, one always obtains < 0. 
Hence the roots are always real and three real unequal roots can be computed as: 

'rh = 2V^cos f - j - -j-, \h - 2y^cos ( j - y, 



= 2-\/— Sicos 



e + 47r 



£1 
3 



(69) 



One finds that for all 1 ^ {^,1} ^ 2, becomes negative. It is observed that {^h,^rh} > 1 for most values of the 
astrophysically tuned {£, 7}. However, it is also found that does not allow steady physical flows to pass through it; 
either u, or Cs, or both, becomes superluminal before the flow reaches the actual event horizon, or the Mach number 
profile shows intrinsic fluctuations for r < r/j. This information is obtained by numerically integrating the complete 
flow profile passing through ^r^. Hence it turns out that one needs to concentrate only on for realistic astrophysical 
black hole accretion. Both large £ and large 7 enhance the thermal energy of the flow so that the accreting fluid 
acquires the large radial velocity to overcome as only in the close vicinity of the black hole . Hence anti-correlates 
with {£,7}. 

The critical properties and stability of such acoustic horizons has recently been studied using a dynamical systems 
approach, see Mandal, Ray & Das 2007 for further details. 

To obtain (du/dr) and [dcg/dr) on the acoustic horizon, L' Hospital's rule is applied to (|63p to have: 



^12 - ^'123, ( 
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$12 *123 
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(70) 



where 



*12 = -$2/2$!, $123 = V^'i - 4$i$3/2$i, $1 = 



(4r,,-3)^ 



Qru {rh - 1) 



$2 = 



4r,, - 3 



[4r,(7-l)-(37-2)] 



rl (7 - 1)^ - rir (107^ - 197 + 9) + (672 - II7 + 3) 



$4 



2 {2rh - 1) - 7 (4rh - 3) 
4 {rh - 1) 



(71) 



C. Analogue Temperature 

For spherically symmetric general relativistic flow onto Schwarzschild black holes, one can evaluate the exact value 
of the Killing fields and Killing vectors to calculate the surface gravity for that geometry. The analogue Hawking 
temperature for such geometry comes out to be (Das 2004a) 



Ta 



H 



AttkbGMbh 



r^^ {rh - 0.75) 

{rn - 1) 



3/2 



dr 



{Cs - U) 



(72) 



where the values of r^, {du/dr)h and {dcs/dr)h are obtained using the system of units and scaling used in this article. 

It is evident from (|72|) that the exact value of Tah can be analytically calculated from the results obtained in 
the previous section. While provides the location of the acoustic horizon (r/j), the value of 1-^ (c — u)| _ is 
obtained from (|70ll7ip as a function of £ and 7, both of which are real, physical, measurable quantities. Note again, 
that, since rh and other quantities appearing in (|72p are analytically calculated as a function of {£,7}, (j72p provides 
an exact analytical value of the general relativistic analogue Hawking temperature for all possible solutions of an 
spherically accreting astrophysical black hole system, something which has never been done in the literature before. If 
\/4,rh — 3(1/2— l/$4)($i2 — $123) > 1, one always obtains {dcg/dr < du/dr)h from ((70)) . which indicates the presence 
of the acoustic white holes at rh- This inequality holds good for certain astrophysically relevant range of {£, 7}, thus 
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acoustic white hole solutions are obtained for general relativistic, spherically symmetric black hole accretion, see Das 
2004a for further detail. 

For a particular value of {^,7}, one can define the quantity r to be the ratio of Tah and Th as: 

r=^. (73) 

It turns out that r is independent of the mass of the black hole. Thus, by computing the value of r, we can compare 
the properties of the acoustic versus event horizon of an accreting black hole of any mass, starting from the primordial 
black hole to the super massive black holes at the dynamical centre of the galaxies. 

For general relativistic spherical accretion, one finds that for certain range of {£,7}, T^u exceeds (i.e., r > 1) the 
value of Tff , hence the analogue Hawking temperature can be larger than the actual Hawking temperature, see Das 
2004a for further details. 



XIII. MULTI-TRANSONIC, RELATIVISTIC ACCRETION DISC AS ANALOGUE GRAVITY MODEL 

A. The Stress Energy Tensor and Flow Dynamics 

To provide a generalized description of axisymmetric fluid flow in strong gravity, one needs to solve the equations of 
motion for the fluid and the Einstein equations. The problem may be made tractable by assuming the accretion to be 
non-self gravitating so that the fluid dynamics may be dealt in a metric without back-reactions. To describe the flow, 
we use the Boyer-Lindquist co-ordinate (Boyer & Lindquist 1967), and an azimuthally Lorentz boosted orthonormal 
tetrad basis co-rotating with the accreting fluid. We define A to be the specific angular momentum of the flow. 
Since we are not interested in non-axisymmetric disc structure, we neglect any gravo-magneto- viscous non-alignment 
between A and black hole spin angular momentum. We consider the flow to be non-self-gravitating to exclude any 
back reaction on the metric. For this section, the gravitational radius rg is taken to be GMbh/(?- 

The most general form of the energy momentum tensor for the compressible hydromagnetic astrophysical fluid 
(with a frozen in magnetic held) vulnerable to the shear, bulk viscosity and generalized energy exchange, may be 
expressed as (Novikov & Thorne 1973): 

c.^,u = cjAj^ + cj^^ (74) 

where and are the fluid (matter) part and the Maxwellian (electromagnetic) part of the energy momentum 
tensor. 9^ and ?s'^ may be expressed as: 

= pv^'v'' + {p- ^e) /i^'' - 2ria^"' + q^'v" + i;^q^ = ^ (B^v'*?;^ + B'^h"'' - 2^"^") (75) 

In the above expression, pv'^v" is the total mass energy density excluding the frozen-in magnetic held mass energy 
density as measured in the local rest frame of the baryons (local orthonormal frame, hereafter LRF, in which there is 
no net baryon flux in any direction), ph'^^ is the isotropic pressure for incompressible gas (had it been the case that 
9 would be zero). <, and ry are the co-efficient of bulk viscosity and of dynamic viscosity, respectively. Hence —(,0h^'^ 
and —2r]a^" are the isotropic viscous stress and the viscous shear stress, respectively, ofv'^ + v^(f is the energy and 
momentum flux, respectively, in LRF of the baryons. In the expression for ^g"^, B^/Stt in the flrst term represents 
the energy density, in the second term represents the magnetic pressure orthogonal to the magnetic field lines, and in 
third term magnetic tension along the; field lines (all terms expressed in LRF), respectively. 

Here, the electromagnetic field is described by the field tensor and it's dual J^*!^" (obtained from Ti^" using 
Levi-Civita 'fiipping' tensor e'^""^) satisfying the Maxwell equations through the vanishing of the four-divergence of 

, A complete description of flow behaviour could be obtained by taking the co- variant derivative of SJ'*'' and 
to obtain the energy momentum conservation equations and the conservation of baryonic mass. 

However, at this stage, the complete solution remains analytically untenable unless we are forced to adopt a 
number of simplified approximations. We would like to study the inviscid accretion of hydrodynamic fluid. Hence 
may be described by the standard form of the energy momentum (stress-energy) tensor of a perfect perfect fiuid: 

'^'"^ = {€+p)v^,v„+pg^,^, or, T = (e+p)v(g)v-Fpg (76) 

Our calculation will thus be focused on the stationary axisymmetric solution of the energy momentum and baryon 
number conservation equations 



S'^^. = 0; {pv^.^ = 0, 



(77) 
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Specifying the metric to be stationary and axially symmetric, the two generators ^'^ = {d/dt)^ and (/>'^ = {d/dcf))^ of 
the temporal and axial isometry, respectively, are Killing vectors. 

We consider the flow to be 'advective', i.e. to possess considerable radial three- velocity. The above-mentioned 
advective velocity, which we hereafter denote by u and consider it to be confined on the equatorial plane, is essentially 
the three-velocity component perpendicular to the set of hypersurfaces {St,} defined by — const, where v is the 
magnitude of the 3-velocity. Each is timelike since its normal rj^ oc d^v^ is spacelike and may be normalized as 
?7''77^ = 1. 

We then define the specific angular momentum A and the angular velocity U, as 



A 



gt4> + Agtt 
+ ^9t<p 



The metric on the equatorial plane is given by (Novikov & Thorne 1973) 



A 9 
dt'^ + ^ (d</« - ujdt) + —dr'^ 
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where A = — 2r 



A = 



hole spin. The normalization condition 



2ra^, and oj = 2ar/A, a being the Kerr parameter related to the black- 
Vn = —1, together with the expressions for A and H. in ([75)) . provides the 



relationship between the advective velocity u and the temporal component of the four velocity 
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(1 - 7/2) {^2 _ ^XarA + A2r2 {Aa^ - r2A)} 
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(80) 



In order to solve (|77p . we need to specify a realistic equation of state. In this work, we concentrate on polytropic 
accretion. However, polytropic accretion is not the only choice to describe the general relativistic axisymmetric black- 
hole accretion. Equations of state other than the adiabatic one, such as the isothermal equation (Yang & Kafatos 
1995) or the two-temperature plasma (Manmoto 2000), have also been used to study the black-hole accretion flow. 

Like spherical accretion, here also we assume the dynamical in-fall time scale to be short compared with any 
dissipation time scale during the accretion process. We have performed the calculations for a sufficiently large range 
of 7 and we believe that all astrophysically relevant polytropic indices are covered in our work. 



B. Disc Geometry and the Conservation Equations 

We assume that the disc has a radius-dependent local thickness H{r), and its central plane coincides with the equatorial 
plane of the black hole. It is a standard practice in accretion disc theory (Matsumoto et. al. 1984; Paczy'nski 1987; 
Abramowicz, Czerny, Lasota & Szuszkiewicz 1988; Chen & Taam 1993; Kafatos & Yang 1994; Artemova, Bjornsson & 
Novikov 1996; Narayan, Kato & Honma 1997; Wiita 1999; Hawley & Krolik 2001; Armitage, Reynolds & Chiang 2001) 
to use the vertically integrated model in describing the black-hole accretion discs where the equations of motion apply 
to the equatorial plane of the black hole, assuming the flow to be in hydrostatic equilibrium in the transverse direction. 
The assumption of hydrostatic equilibrium is justifled for a thin flow because for such flows, the infall time scale is 
expected to exceed the local sound crossing time scale in the direction transverse to the flow. We follow the same 
procedure here. The thermodynamic flow variables are averaged over the disc height, i.e., a thermodynamic quantity 
y used in our model is vertically integrated over the disc height and averaged as y = {r)^{ydh) / {r)QH{r). 



In flgure 2, we schematically represent the above mentioned modelling. The yellow circular patch with BH written 
inside represents the black hole and the pink dashed boundary mimics the event horizon. The wedge shaped dark 
green lines represents the envelop of the accretion disc. The light green line centrally flanked by the two dark green 
disk boundaries, is the equatorial plane, on which all of the dynamical quantities (e.g., the advective velocity u) are 
assumed to be confined. Any thermodynamic quantity (e.g., the flow density) is averaged over the local disc height 
h(r) as shown in the figure. 

We follow Abramowicz, Lanza & Percival (1997) to derive an expression for the disc height H{r) in our fiow 
geometry since the relevant equations in Abramowicz, Lanza & Percival (1997) are non-singular on the horizon and 
can accommodate both the axial and a quasi-spherical flow geometry. In the Newtonian framework, the disc height 
in vertical equilibrium is obtained from the z component of the non-relativistic Euler equation where all the terms 
involving velocities and the higher powers of (z/r) are neglected. In the case of a general relativistic disc, the vertical 
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FIG. 2: Height averaged thermodynamic quantities for disc accretion. 



pressure gradient in the comoving frame is compensated by the tidal gravitational field. We then obtain the disc 
height 



Hir) = 



7+1 



{j~{l + cf}}{X^vf-a^(vt-l)} 



(81) 



which, by making use of (j80p . may be be expressed in terms of the advective velocity u. 

The temporal component of the energy momentum tensor conservation equation leads to the constancy along 
each streamline of the flow specific energy £ {£ = hvt), and hence from (|80p and ([55)1 it follows that: 
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The rest- mass accretion rate M is obtained by integrating the relativistic continuity equation (|77p . One finds 

1 



(82) 
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Here, we adopt the sign convention that a positive u corresponds to accretion. The entropy accretion rate can be 
expressed as: 
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One can solve the conservation equations for £, M and to obtain the complete accretion profile. 



C. Transonicity 

The gradient of the acoustic velocity can be computed by differentiating ([84)1 and can be obtained as: 
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The dynamical velocity gradient can then be calculated by differentiating (ISS)) with the help of ((85|) as: 



du 
dr 



2c? 


V - 1 ^ 2 wtCTx' 


X 


(7+1) 


A r 4^^ 


2 


u 


2c? 


1 




(1-m2) (7+1)(1-u2)u 


2V' 



(86) 



where 
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The critical point conditions can be obtained as: 
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For any value of [£ , A, 7, a], substitution of the values of u|(r=rc) s-'^d Cs|r=rc in terms of rc in the expression for £ (j82p . 

provides a polynomial in rc, the solution of which determines the location of the critical point(s) ra- 
il is obvious from that, unlike relativistic spherical accretion, )7^Csj.=ro; ^^id hence the Mach number 

at the critical point is not equal to unity in general. This phenomena can more explicitly be demonstrated for = 0, 

i.e., for relativistic disc accretion in the Schwarzschild metric. 

For Schwarzschild black hole, one can calculate the Mach number of the flow at the critical point as j228j | (Das, 

Bflic & Dasgupta 2006) 
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Clearly, Mj. is generally not equal to unity, and for 7 > 1, is always less than one. 

Hence we distinguish a sonic point from a critical point. In the literature on transonic black-hole accretion discs, 
the concepts of critical and sonic points are often made synonymous by defining an 'effective' sound speed leading to 
the 'effective' Mach number (for further details, see, eg. Matsumoto et. al. 1984, Chakrabarti 1989). Such definitions 
were proposed as effects of a specific disc geometry. We, however, prefer to maintain the usual definition of the Mach 
number for two reasons. 

First, in the existing literature on transonic disc accretion, the Mach number at the critical point turns out to be 
a function of 7 only, and hence Mc remains constant if 7 is constant. For example, using the Paczyhski and Wiita 
(1980) pseudo-Schwarzschild potential to describe the adiabatic accretion phenomena leads to (see section 16.1.1 for 
the derivation and for further details) 



Mr 



7+1 



(91) 



The above expression does not depend on the location of the critical point and depends only on the value of the 
adiabatic index chosen to describe the flow. Note that for isothermal accretion 7 = 1 and hence the sonic points and 
the critical points are equivalent (since Mc = 1), see (jl69|) in section 16.1.2 for further details. 

However, the quantity Mc in Eq. ([55]) as well as in ((55|) is clearly a function of r^ and hence, generally, it takes 
different values for different rc for transonic accretion. The difference between the radii of the critical point and the 
sonic point may be quite significant. One defines the radial difference of the critical and the sonic point (where the 
Mach number is exactly equal to unity) as 



Ar^ = \rs - rc\. 



(92) 
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The quantity Ar^ may be a complicated function of [£", A,7,a], the form of which can not be expressed analytically. 
The radius in Eq. ([M]) is the radius of the sonic point corresponding to the same [£, A, 7, a]for which the radius of 
the critical point rj, is evaluated. Note, however, that since is calculated by integrating the flow from Tc, Ar* is 
defined only for saddle-type critical points (see subsequent paragraphs for further detail). This is because, a physically 
acceptable transonic solution can be constructed only through a saddle- type critical point. One can then show that 
Ar^ can be as large as 10^ rg or even more (for further details, see Das, Bilic & Dasgupta 2006). 

The second and perhaps the more important reason for keeping Vc and Vs distinct is the following. In addition 
to studying the dynamics of general relativistic transonic black-hole accretion, we are also interested in studying the 
analogue Hawking effects for such accretion flow. We need to identify the location of the acoustic horizon as a radial 
distance at which the Mach equals to one, hence, a sonic point, and not a critical pomi will be of our particular interest. 
To this end, we first calculate the critical point Tc for a particular [f , A, 7, aJfoUowing the procedure discussed above, 
and then we compute the location of the sonic point (the radial distance where the Mach number exactly equals to 
unity) by integrating the flow equations starting from the critical points. The dynamical and the acoustic velocity, 
as well as their space derivatives, at the sonic point, are then evaluated. The details of this procedure for the 
Schwarzschild metric are provided in Das, Bilic & Dasgupta 2006. 

Furthermore, the definition of the acoustic metric in terms of the sound speed does not seem to be mathemat- 
ically consistent with the idea of an 'effective' sound speed, irrespective of whether one deals with the Newtonian, 
post Newtonian, or a relativistic description of the accretion disc. Hence, we do not adopt the idea of identifying 
critical a point with a sonic point. However, for saddle- type critical points, and should always have one-to-one 
correspondence, in the sense that every critical point that allows a steady solution to pass through it is accompanied 
by a sonic point, generally at a different radial distance r. 

It is worth emphasizing that the distinction between critical and sonic points is a direct manifestation of the 
non-trivial functional dependence of the disc thickness on the fluid velocity, the sound speed and the radial distance. 
In the simplest idealized case when the disc thickness is assumed to be constant, one would expect no distinction 
between critical and sonic points. In this case, as has been demonstrated for a thin disc accretion onto the Kerr black 
hole (Abraham, Bilic & Das 2006), the quantity Ar^ vanishes identically for any astrophysically relevant value of 
[f , A,7,a]. Hereafter, we will use to denote the sonic point rg, since a sonic point is actually the location of the 
acoustic horizon. 
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FIG. 3: Parameter space for general relativistic multi-transonic accretion and wind in Kerr geometry, see text for detail. This 
figure is reproduced from Goswami, Khan, Ray & Das 2007. 
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D. Multi-transonic Behaviour 



Unlike spherical accretion, one finds three (at most) critical points for relativistic disc accretion for some values of 
[£, A, 7, a]. In figure 3, we classify the [£, A] parameter space, for a fixed value of adiabatic index (7 = 4/3) and the 
Kerr parameter (a = 0.3), to show the formation of various kind of critical points. The regions marked by O and I 
correspond to the formation of a single critical point, and hence the mono-transonic disc accretion is produced for 
such region. In the region marked by I, the critical points are called 'inner type' critical points since these points 
are quite close to the event horizon, approximately in the range 2 < r™<lO. In the region marked by O, the critical 
points are called 'outer type' critical points, because these points are located considerably far away from the black 
hole. Depending on the value of [£, A, 7, a], an outer critical point may be as far as lO^r^, or more. 

The outer type critical points for the mono-transonic region are formed, as is obvious from the figure, for weakly 
rotating flow. For low angular momentum, accretion flow contains less amount of rotational energy, thus most of 
the kinetic energy in utilized to increase the radial dynamical velocity m at a faster rate, leading to a higher value of 
cPu/dr'^ . Under such circumstances, the dynamical velocity u becomes large enough to overcome the acoustic velocity 
Cs at a larger radial distance from the event horizon, leading to the generation of supersonic flow at a large value of 
r, which results the formation of the sonic point (and hence the corresponding critical point) far away from the black 
hole event horizon. On the contrary, the inner type critical points are formed, as is observed from the figure, for 
strongly rotating flow in general. Owing to the fact that such flow would posses a large amount of rotational energy, 
only a small fraction of the total specific energy of the flow will be spent to increase the radial dynamical velocity 
u. Hence for such flow, u can overcome Cs only at a very small distance (very close to the event horizon) where the 
intensity of the gravitational field becomes enormously large, producing a very high value of the linear kinetic energy 
of the fiow (high u), over shedding the contribution to the total specific energy from all other sources. However, from 
the figure it is also observed that the inner type sonic points are formed also for moderately low values of the angular 
momentum as well (especially in the region close to the vertex of the wedge shaped zone marked by W). For such 
regions, the total conserved specific energy is quite high. In the asymptotic limit, the expression for the total specific 
energy is governed by the Newtonian mechanics and one can have: 



^^(y + Witational (93) 

/ linear V ' / thermal V / rotational 



where $ is the gravitational potential energy in the asymptotic limit, see section 16.1.1 for further detail. From ((93|) 
it is obvious that at a considerably large distance from the black hole, the contribution to the total energy of the flow 
comes mainly (rather entirely) from the thermal energy. A high value of E (flow energy in excess to its rest mass 
energy) corresponds to a 'hot' flow starting from infinity. Hence the acoustic velocity corresponding to the 'hot' flow 
obeying such outer boundary condition would be quite large. For such accretion, flow has to travel a large distance 
subsonically and can acquire a supersonic dynamical velocity u only at a very close proximity to the event horizon, 
where the gravitational pull would be enormously strong. 

The [£, A] corresponding to the wedge shaped regions marked by A and W produces three critical points, among 
which the largest and the smallest values correspond to the X type (saddle type), the outer r°"* and the inner r™, 
critical points respectively. The O type (centre type) middle critical point, r™*"*, which is unphysical in the sense that 
no steady transonic solution passes through it, lies in between r™ and r°"*. The following discussion provides the 
methodology for finding out the nature (whether saddle/centre type) in brief, see Goswami, Khan, Ray & Das 2007 
for further detail. 
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with the primes representing full derivatives with respect to r, and f is an arbitrary mathematical parameter. Here, 
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The critical conditions are obtained with the simultaneous vanishing of the right hand side, and the coefficient of 
d{u'^)/dr in the left hand side in ([M]) . This will provide 
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as the two critical point conditions. Some simple algebraic manipulations will show that 
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following which c^jr^i-^ can be rendered as a function of Tc only, 
be fixed in terms of the constants of motion like E, 7, A and a. 
necessary to study their nature in their phase portrait of vP' versus r. To that end one applies a perturbation about 
the fixed point values, going as, 
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in the parameterized set of autonomous first-order differential equations. 
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7 + 1 ^ 

with f being an arbitrary parameter. In the two equations above Jc^ can be closed in terms of Sv?' and 5r with the 
help of ([55]) . Having done so, one could then make use of solutions of the form, 5r ~ exp(r2r) and Sv?' ~ exp(r2r), 
from which, f2 would give the eigenvalues — growth rates of 5u^ and 5r in f space — of the stability matrix implied 
by p00m01|) . Detailed calculations will show the eigenvalues to be 

f^' = |/3'4xi+Ci6L=,^ (102) 

where pi^ = and Xi^^i and ^2 can be expressed as polynomials of Tc (see Goswami, Khan, Ray & Das 2007 for 
the explicit form of the polynomial), hence CP' can be evaluated for any [£, A, 7, ajonce the value of the corresponding 
critical point Tc is known. The structure of (|102p immediately shows that the only admissible critical points in the 
conserved Kerr system will be either saddle points or centre type points. For a saddle point, > 0, while for a 
centre- type point, < 0. 

For multi-transonic flow characterized by a specific set of [£, A, 7, a], one can obtain the value of Cl^ to be positive 



for 



and ' 



showing that those critical points are of saddle type in nature, comes out to be negative for 



confirming that the middle sonic point is of centre type and hence no transonic solution passes through it. One can 
also confirm that all mono-transonic flow (flow with a single critical point characterized by [£, A] used from the green 
tinted region, either I or O) corresponds to saddle type critical point. 

However, there is a distinct difference between the multi-transonic flow characterized by [£, A] taken from the 
region marked by A, and the region marked by W. For region marked by A, the entropy accretion rate S for flows 
passing through the inner critical point is greater than that of the outer critical point 

(103) 



while for the region marked by W, the following relation holds 

S (rT) < E (r° 



(104) 



The above two relations show that [£, A] region marked by A represents multi-transonic accretion, while [£, A] G [£, A]-^ 
corresponds to the mono-transonic accretion but multi-transonic wind. More details about such classification will be 
discussed in the following paragraphs. 

There are other regions for [£, A] space for which either no critical points are formed, or two critical points are 
formed. These regions are not shown in the figure. However, none of these regions is of our interest. If no critical point 
is found, it is obvious that transonic accretion does not form for those set of , A]. For two critical point region, one 
of the critical points are always of 'O' type, since according to the standard dynamical systems theory two successive 
critical points can not be of same type (both saddle, or both centre). Hence the solution which passes through the 
saddle type critical point would encompass the centre type critical point by forming a loop (see, e.g.. Das, Bilic & 
Dasgupta 2006 for such loop formation in Schwarzschild metric) like structure and hence such solution would not be 
physically acceptable since that solution will form a closed loop and will not connect infinity to the event horizon. 
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E. Multi-transonic Flow Topology and Shock Formation 



To obtain the dynamical velocity gradient at the critical point, we apply I'Hospital's rule on ((86|) . After some algebraic 
manipulations, the following quadratic equation is formed, which can be solved to obtain {du/ dr)]^^^_^ ^ (see Barai, 
Das & Wiita 2004 for further details): 



where the coefficients are: 
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Note that all the above quantities are evaluated at the critical point. 
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Hence we compute the critical advective velocity gradient as 

where the '+' sign corresponds to the accretion solution and the '-' sign corresponds to the wind solution, see the 
following discussion for further details. Similarly, the space gradient of the acoustic velocity dcs/dr and its value at 
the critical point can also be calculated. 

The flow topology characterized by [£", A] corresponding to the I or O region (green tinted) is easy to obtain since 
the flow passes through only one saddle type critical point. Some of such representative topologies will be found in 
Das, Bilic & Dasgupta 2006 for the Schwarzschild metric. For Kerr metric, the flow profile would exactly be the 
same, only the numerical values for the critical/sonic point would be different, leading to the different values of u, 
Cs and other corresponding accretion parameters at the same radial distance r. In this section we concentrate on 
multi-transonic flow topology, i.e., flow topology for [£, A] G , A]^ or [£, A] G [5, A].^. Figure 4 represents one such 
topology. While the local radial flow Mach number has been plotted along the Y axis, the distance from the event 
horizon (scaled in the unit of GMbh / <?) in logarithmic unit has been plotted along the X axis. 

The solid red line marked by A corresponds to the transonic accretion passing through the outer critical point 
r°"*. Using a specific set of [£,A,7,a]as shown in the figure, one first solve the equation for £ at the critical point 
(using the critical point conditions as expressed in (jSSp to find out the corresponding three critical points, saddle 
type r™ (4.279 r^), centre type r"''' (14.97 r^) and saddle type (3315.01 Vg). The critical value of the advective 
velocity gradient at r°"* is then calculated using (|105m07|) . Such u\^^^_^^^,Cs\^^^^^^ and du/ dr\^^^^^^ serve as the initial 
value condition for performing the numerical integration of the advective velocity gradient (|86p using the fourth-order 
Runge-Kutta method. Such integration provides the outer sonic point r°"* (r°"* < r°"*), the local advective velocity, 
the polytropic sound speed, the Mach number, the fluid density, the disc height, the bulk temperature of the flow, 
and any other relevant dynamical and thermodynamic quantity characterizing the flow. 

The dotted green line marked by W is obtained for the value of du/dr\^^_^ ^ corresponds to the '-' sign in (|107p . 
Such a solution is called the corresponding 'wind' solution. The accretion and the wind solution intersects at the 
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FIG. 4: Solution topology for multi-transonic accretion in Kerr geometry for a specific set of [£■, A,7,a]as shown in the figure. 
See text for detail. 



critical point (here, at r°"*). This wind branch is just a mathematical counterpart of the accretion solution (velocity 
reversal symmetry of accretion) , owing to the presence of the quadratic term of the dynamical velocity in the equation 
governing the energy momentum conservation. The term 'wind solution' has a historical origin. The solar wind 
solution first introduced by Parker (1965) has the same topology profile as that of the wind solution obtained in 
classical Newtonian Bondi accretion (Bondi 1952). Hence the name 'wind solution' has been adopted in a more 
general sense. The wind solution thus represents a hypothetical process, in which, instead of starting from infinity 
and heading towards the black hole, the flow generated near the black-hole event horizon would fly away from the 
black hole towards infinity. 

The dashed blue line marked by 'a' and the dotted magenta line marked by 'w' are the respective accretion 
and the wind solutions passing through the inner critical point (the intersection of the accretion and the wind 
branch is the location of r™). Such accretion and wind profile are obtained following exactly the same procedure 
as has been used to draw the accretion and wind topologies (red and green lines) passing through the outer critical 
point. Note, however, that the accretion solution through r^" folds back onto the wind solution and the accretion-wind 
closed loop encompasses the middle sonic point, location of which is represented by M in the figure. One should note 
that an 'acceptable' physical transonic solution m ust b e globally consistent, i.e. it must connect the radial infinity 
r— !-oo with the black-hole event horizon r = 2rg. [229( . Hence, for multi-transonic accretion, there is no individual 
existence of physically acceptable accretion/ wind solution passing through the inner critical (sonic) point, although 
such solution can be 'clubbed' with the accretion solution passing through r°"* through shock formation, see the 
following discussions for further details. 

The set , A]^ (or more generally [f , A, 7, a]^) thus produces doubly degenerate accretion/ wind solutions. Such 
two fold degeneracy may be removed by the entropy considerations since the entropy accretion rates S(r™) and S(r°"*) 
are generally not equal. For any [£, A, 7, a] £ [£, A, 7, a]^ we find that the entropy accretion rate S evaluated for the 
complete accretion solution passing through the outer critical point is less than that of the rate evaluated for the 
incomplete accretion/ wind solution passing through the inner critical point. Since the quantity S is a measure of 
the specific entropy density of the flow, the solution passing through r°"* will naturally tend to make a transition 
to its higher entropy counterpart, i.e. the incomplete accretion solution passing through r^". Hence, if there existed 
a mechanism for the accretion solution passing through the outer critical point (solid red line marked with A) to 
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increase its entropy accretion rate by an amount 

AS = S(C)-5(rr), (108) 

there would be a transition to the incomplete accretion solution (dashed blue line marked with 'a') passing through 
the inner critical point. Such a transition would take place at a radial distance somewhere between the radius of the 
inner sonic point and the radius of the accretion/ wind turning point (75.7 Vg) marked by G in the figure. In this 
way one would obtain a combined accretion solution connecting r— s-oo with r = 2 (the event horizon) which includes 
a part of the accretion solution passing through the inner critical, and hence the inner sonic point. One finds that 
for some specific values of [£, A,7,a]^, a standing Rankine-Hugoniot shock may accomplish this task. A supersonic 
accretion through the outer sonic point r°"' (which in obtained by integrating the flow starting from the outer critical 
point r°"*) can generate entropy through such a shock formation and can join the flow passing through the inner 
sonic point r™ (which in obtained by integrating the flow starting from the outer critical point r™). Below we will 
carry on a detail discussion on such shock formation. 

In this article, the basic equations governing the flow are the energy and baryon number conservation equations 
which contain no dissipative terms and the flow is assumed to be inviscid. Hence, the shock which may be produced 
in this way can only be of Rankine-Hugoniot type which conserves energy. The shock thickness must be very small 
in this case, otherwise non-dissipative flows may radiate energy through the upper and the lower boundaries because 
of the presence of strong temperature gradient in between the inner and outer boundaries of the shock thickness. In 
the presence of a shock the flow may have the following profile. A subsonic flow starting from infinity first becomes 
supersonic after crossing the outer sonic point and somewhere in between the outer sonic point and the inner sonic 
point the shock transition takes place and forces the solution to jump onto the corresponding subsonic branch. The 
hot and dense post-shock subsonic flow produced in this way becomes supersonic again after crossing the inner sonic 
point and ultimately dives supersonically into the black hole. A flow heading towar ds a neutron star can have the 
liberty of undergoing another shock transition after it crosses the inner sonic point '230], because the hard surface 
boundary condition of a neutron star by no means prevents the flow from hitting the star surface subsonically. 

For the complete general relativistic accretion flow discussed in this article, the energy momentum tensor 5'^", 
the four- velocity v^, and the speed of sound may have discontinuities at a hypersurface E with its normal 77^. Using 
the energy momentum conservation and the continuity equation, one has 

[[p«^]]77^ = 0,[[3n]'?. = 0. (109) 

For a perfect fluid, one can thus formulate the relativistic Rankine-Hugoniot conditions as 

[[puTu]] = 0, (110) 



[[%^r]]^[[iP + e)vtuT^]]^0, (111) 

[[^^..rvl] ^[[{P + ^yrl + p] ] = 0, (112) 

where r,j — l/Vl — is the Lorentz factor. The first two conditions ()110|1 and (llll|) are trivially satisfied owing to 
the constancy of the specific energy and mass accretion rate. The constancy of mass accretion yields 

_[ \ 1 J W-i-c^; 713^ ^ > 

The third Rankine-Hugoniot condition (jll2p may now be written as 

[ \ 1 ) \l-l-cl) \ cl{l-u^) / 

Simultaneous solution of Eqs. pisp and (|114p yields the 'shock invariant' quantity 

Su = cP^ (7 - 1 - c^)^ (1 - u'Y'^ [Wl - a" {vt - 1)]"^ 

which changes continuously across the shock surface. We also define the shock strength Si and the entropy enhancement 
Q as the ratio of the pre-shock to post-shock Mach numbers {Si = M^/M+), and as the ratio of the post-shock to 
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pre-shock entropy accretion rates {Q = S+/S_) of the flow, respectively. Hence, Q = S(r*")/Sr°'"* for accretion and 
Q = S(r°"*)/Sr™ for wind, respectively. 

The shock location in multi-transonic accretion is found in the following way. Consider the multi-transonic flow 
topology as depicted in the Fig. 4. Integrating along the solution passing through the outer critical point, we calculate 
the shock invariant Sh in addition to u, Cs and M . We also calculate Sh while integrating along the solution passing 
through the inner critical point, starting from the inner sonic point up to the point of inflexion G. We then determine 
the radial distance rsh , where the numerical values of Sh , obtained by integrating the two different sectors described 
above, are equal. Generally, for any value of [£, A, 7, ajallowing shock formation, one finds two shock locations marked 
by SI (the 'outer' shock, formed at 65.31rg - between the outer and the middle sonic points) and S2 (the 'inner' 
shock, formed at 6.31 - between the inner and the middle sonic points) in the figure. According to a standard local 
stability analysis (Yang & Kafatos 1995), for a multi-transonic accretion, one can show that only the shock formed 
between the middle and the outer sonic point is stable. The shock strength is different for the inner and for the outer 
shock. For the stable (outer) shock, the shock strength for the case shown in the figure is 5.586, hence it is a strong 
shock. Therefore, in the multi-transonic accretion with the topology shown in Fig. 4., the shock at SI is stable and 
that at S2 is unstable. Hereafter, whenever we mention the shock location, we refer to the stable shock location only. 
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FIG. 5: Pre- and post-shock disc geometry with thermally driven optically thick halo. See text for further detail. 



F. Disc Geometry and Shock Generated Outflow 



As a consequence of the shock formation in an accretion flow, the post-shock flow temperature will also increase 
abruptly. The post- to pre-shock temperature ratio T^jT^ is, in general, a sensitive function of [£, A,7,a]. In Fig. 
5., we present the disc structure obtained by solving ((8T|) for the combined shocked accretion flow. The point B 
represents the black-hole event horizon. The pre- and post-shock regions of the disc are clearly distinguished in the 
figure and show that the post-shock disc puffs up significantly. The pre-shock supersonic disc is shown by the red 
lines. The post-shock subsonic part of the disc is shown by dotted blue lines and the post-shock supersonic part 
(very close to the event horizon since r™ = 4.279rg) is shown by dotted magenta lines (not well resolved in the figure 
though). The bulk flow temperature will be increased in the post-shock region. Such an increased disc temperature 
may lead to a disc evaporation resulting in the formation of an optically thick halo, which are schematically shown 
using yellow coloured elliptic structures. Besides, a strong temperature enhancement may lead to the formation 
of thermally driven outfiows. The generation of centrifugally driven and thermally driven outflows from black-hole 
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accretion discs has been discussed in the post-Newtonian framework (Das & Chakrabarti 1999; Das, Rao & Vadawale 
2003) The post-Newtonian approach may be extended to general relativity using the formalism presented here. 

Owing to the very high radial component of the infall velocity of accreting material close to the black hole, the 
viscous time scale is much larger than the infall time scale. Hence, in the vicinity of the black hole, a rotating inflow 
entering the black hole will have an almost constant specific angular momentum for any moderate viscous stress. This 
angular momentum yields a very strong centrifugal force which increases much faster than the gravitational force. 
These two forces become comparable at some specific radial distance. At that point the matter starts piling up and 
produces a boundary layer supported by the centrifugal pressure, which may break the inflow to produce the shock. 
This actually happens not quite at the point where the gravitational and centrifugal forces become equal but slightly 
farther out owing to the thermal pressure. Still closer to the black hole, gravity inevitably wins and matter enters the 
horizon supersonically after passing through a sonic point. The formation of such a layer may be attributed to the 
shock formation in accreting fluid. The post-shock flow becomes hotter and denser, and for all practical purposes, 
behaves as the stellar atmosphere as far as the formation of outflows is concerned. A part of the hot and dense 
shock-compressed in-flowing material is then 'squirted' as an outflow from the post-shock region. Subsonic outflows 
originating from the puffed up hotter post-shock accretion disc (as shown in the figure) pass through the outflow sonic 
points and reach large distances as in a wind solution. 

The generation of such shock-driven outflows is a reasonable assumption. A calculation describing the change of 
linear momentum of the accreting material in the direction perpendicular to the plane of the disc is beyond the scope 
of the disc model described in this article because the explicit variation of dynamical variables along the Z axis (axis 
perpendicular to the equatorial plane of the disc) cannot be treated analytically. The enormous post-shock thermal 
pressure is capable of providing a substantial amount of 'hard push' to the accreting material against the gravitational 
attraction of the black hole. This 'thermal kick' plays an important role in re-distributing the linear momentum of 
the inflow and generates a non-zero component along the Z direction. In other words, the thermal pressure at the 
post-shock region, being anisotropic in nature, may deflect a part of the inflow perpendicular to the equatorial plane 
of the disc. Recent work shows that (Moscibrodzka, Das & Czerny 2006) such shock-outflow model can be applied to 
successfully investigate the origin and dynamics of the strong X-ray flares emanating out from our galactic centre. 

G. Multi-transonic Wind 

The blue coloured wedge shaped region marked by W represents the [f. A, 7, ojzone for which three critical points, 
the inner, the middle and the outer are also found. However, in contrast to , A, 7, a] G , A, 7, a]^, the set 
[f, A,7, a] G [£,X,^,a]-^ yields solutions for which S(r™) is less than S(r°"*). Besides, the topological flow profile 
of these solutions is different. Here the closed loop-like structure is formed through the outer critical point. One such 
solution topology is presented in Fig. 6 for a specific set of [£, A, 7, a] as shown in the figure. The same colour-scheme 
which has been used to denote various accretion and wind branches (through various critical points) of multi-transonic 
accretion (Fig. 4.), has been used here as well to manifest how the loop formation switches from flow through r™ 
(multi- transonic accretion, Fig. 4.) to flow through r°"* (multi-transonic wind. Fig. 6.). This topology is interpreted 
in the following way. The flow (blue dashed line marked by 'a') passing through the inner critical point (3.456 rg) 
is the complete mono-transonic accretion flow, and the dotted magenta line marked by 'w' is its corresponding wind 
solution. The solutions passing through the outer critical point (3307.318 rg), represents the incomplete accretion 
(solid red line marked by 'A')/wind (dashed green line marked by 'W') solution. However, as S(r™) turns out to 
be less than S(r°"*), the wind solution through r™ can make a shock transition to join its counter wind solution 
passing through r°"*, and thereby increase the entropy accretion rate by the amount AS = S(r°"*) — S(r™). Here 
the numerical values of Sh along the wind solution passing through the inner critical point are compared with the 
numerical values of Sh along the wind solution passing through the outer critical point, and the shock locations SI 
and S2 for the wind are found accordingly. Here also, two theoretical shock locations are obtained, which are shown 
by dot dashed azure vertical lines marked by SI (at 649.41 rg) and S2 (at 6.42 rg), out of which only one is stable. The 
shock strength corresponding to the stable outer shock can be calculated to be 20.24. Hence extremely strong shocks 
are formed for multi-transonic wind in general. A part of the region [f, A,7,a]Yv thus corresponds to mono-transonic 
accretion solutions with multi-transonic wind solutions with a shock. 

Besides 7 = 4/3 and a = 0.3, for which Fig. 3. has been drawn, one can perform a similar classification for 
any astrophysically relevant value of 7 and a as well. Some characteristic features of [£, A] would be changed as 7 is 
being varied. For example, if Srnax is the maximum value of the energy and if A^ax and Amin are the maximum and 
the minimum values of the angular momentum, respectively, for [£, A]^ for a fixed value of 7, then [£ma.x, Amax, Amin] 
anti-correlates with 7. Hence, as the flow makes a transition from its ultra-relativistic to its purely non-relativistic 
limit, the area representing [£, A]^ decreases. 
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H. Dependence of Shock Location on Accretion Parameters 



One finds that the shock location correlates with A. This is obvious because the higher the flow angular momentum, 
the greater the rotational energy content of the flow. As a consequence, the strength of the centrifugal barrier which 
is responsible to break the incoming flow by forming a shock will be higher and the location of such a barrier will be 
farther away from the event horizon. However, the shock location anti-correlates with £ and 7. This means that for 
the same £ and A, in the purely non-relativistic flow the shock will form closer to the black hole compared with the 
ultra-relativistic flow. Besides, we find that the shock strength Si anti-correlates with the shock location r^/j, which 
indicates that the closer to the black hole the shock forms , the higher the strength Si and the entropy enhancement 
ratio 6 are. The ultra-relativistic flows are supposed to produce the strongest shocks. The reason behind this is also 
easy to understand. The closer to the black hole the shock forms, the higher the available gravitational potential 
energy must be released, and the radial advective velocity required to have a more vigorous shock jump will be larger. 
Besides we note that as the flow gradually approaches its purely non-relativistic limit, the shock may form for lower 
and lower angular momentum, which indicates that for purely non-relativistic accretion, the shock formation may 
take place even for a quasi-spherical flow. However, it is important to mention that a shock formation will be allowed 
not for every [£,A,7, ajs [£,A,7,a]^, Equation (|115p will be satisfied only for a specific subset of [£,A,7,a]_^, for 
which a steady, standing shock solution will be found. 
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I. Analogue Temperature 



The surface gravity is defined according to (|33[) . For axisymmctric accretion described in the above sections, one 
can calculate that (Abraham, Bilic & Das 2006; Das, Bilic & Dasgupta 2006) 



r3 + + 2a2 - 2\a 



where B can be defined as 



Since 
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the expression for the analogue temperature can be calculated as 
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C12 {rh,a, A) = + r^a^ + 2rfta^, 
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rl {a^ - Xa^) + 2rla (a - 2A + 1) + Ar,, (a^ - 2Aa + A)] 
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Using (|84II881 fT05lll07p . along with the expression for {dcs/dr) at rc, one can calculate the location of the acoustic 
horizon (the flow sonic point), and the value of c^, du/dr and dcg/dr at the acoustic horizon, by integrating the flow 
from the critical point upto the acoustic horizon (sonic point). Such values can be implemented in the expression for 
Tah in (|119m20p to calculate the analogue temperature. The ratio t = Tah/Th can also be calculated accordingly. 

One can calculate the analogue temperature for the following five different categories of accretion flow all together, 
since we are not interested at this moment to study the analogue effects in wind solutions: 

1. A/Iono-transonic flow passing through the single inner type critical/sonic point. The range of [£, A, 7, ajused to 
obtain the result for this region corresponds to the region of Fig. 3 marked by I. 

2. Mono-transonic flow passing through the single outer type critical/sonic point. The range of [£, A, 7, ajused to 
obtain the result for this region corresponds to the region of Fig. 3 marked by O. 

3. Multi-transonic accretion passing through the mner critical/sonic point. The range of [£, A, 7, ajused to obtain 
the result for this region corresponds to the region of Fig. 3 marked by A. 

4. Multi-transonic accretion passing through the outer critical/sonic point. The range of [£, A, 7, a]used to obtain 
the result for this region corresponds to the region of Fig. 3 marked by A. 

5. Mono-transonic accretion passing through the inner critical/sonic point for the multi-transonic wind zone. The 
range of [£, A, 7, ajused to obtain the result for this region corresponds to the region of Fig. 3 marked by W. 

In this section we would mainly like to concentrate to the study the dependence of Tah on the Kerr parameter a, 
also, we would like to demonstrate that for some values of [£,A,7, a], the analogue temperature may be comparable 
to the actual Hawking temperature. Hence we are interested in the region of [£, A,7,a]for which r can have a value 
as large as possible. We found that large value of r can be obtained only for very high energy flow with large value 
of the adiabatic index. Such an almost purely nonrelativistic hot accretion does not produce multi-transonicity, it 
produce only mono-transonic flow passing through the inner type critical/sonic point. Hence in the figure 7, we show 
the variation of r with a for a specific value of [£, A, 7] (as shown in the figure) for which [f , A, 7] € [5, A, 7]j. However, 
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same t — a figures can be drawn for any [£, A, 7, ajtaking for any of the other four categories of accretion mentioned 
above. 

In figure 7, the ratio of the analogue to the actual Hawking temperature r has been plotted along the Y axis, 
while the black hole spin parameter (the Kerr parameter a) has been plotted along the X axis. It is obvious from 
the figure that there exists a preferred value of the black hole spin parameter for which the acoustic surface gravity 
attains its maximum value. Location of such a peak in the t — a graph, i.e., the preferred value of the Kerr parameter 
which maximizes the surface gravity, sensitively depends on £, A and 7, see Barai & Das 2007 for further details. This 
is an extremely important finding since it manifest the fact that the black hole spin angular momentum does influence 
the analogue gravity effect, and tells how the background (fluid) metric influences the perturbative (acoustic) metric. 
Note that r > 1 is possible to obtain for an extremely large value of £ having the adiabatic index almost equal to its 
purely non-relativistic limit (7 — 5/3). 
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FIG. 7; Variation of the ratio of analogue to the actual Hawking temperature r with the black hole spin angular momentum 
(the Kerr parameter a). 



As mentioned earlier, the discriminant V of the corresponding acoustic metric changes sign whenever the state 
of transonicity of the flow flips from sub- to super-sonicity or vice versa. For multi-transonic shocked accretion flow, 
such state flipping occurs three times, first, from a sub-sonic to a supersonic state at the outer sonic point (the outer 
acoustic horizon r^"*), then from the supersonic to the subsonic state at the shock location through the Rankine- 
Hugoniot shock, and then from the subsonic to the supersonic state again at the inner sonic point (the inner acoustic 
horizon r™). A transition from V < (subsonic flow) to I? > (supersonic flow) produces an acoustic black hole, 
while the reverse transition {V > — >P < 0) produces an acoustic white hole (Bercelo, Liberati, Sonego & Visser 
2004; Abraham, Bilic & Das 2006). It is thus obvious that for multi-transonic accretion encountering a stable shock, 
two acoustic black holes are formed, one at the inner and the other at the outer acoustic horizons (the inner and 
outer sonic points), and an acoustic white hole is produced at the shock. For relativistic accretion disc with constant 
thickness, this has formally been demonstrated (Abraham, Bilic & Das 2006) by computing the value of D for the 
whole span of r ranging from inflnity to the event horizon to demonstrate that the 2? < — >'D > transition indeed 
takes place at r™ and at r°"*, and V > — >I? < transition takes place at the shock location. Similar calculation 
can also be performed for the disc geometry with the specific form of disc height (|8ip used in this work. 
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XIV. BLACK HOLE ACCRETION IN 'MODIFIED' NEWTONIAN POTENTIALS 



Rigorous investigation of the complete general relativistic multi-transonic black hole accretion disc structure is ex- 
tremely complicated. At the same time it is understood that, as relativistic effects play an important role in the 
regions close to the accreting black hole (where most of the gravitational potential energy is released), purely New- 
tonian gravitational potential cannot be a realistic choice to describe transonic black hole accretion in general. To 
compromise between the ease of handling of a Newtonian description of gravity and the realistic situations described 
by complicated general relativistic calculations, a series of 'modified' Newtonian potentials have been introduced to 
describe the general relativistic effects that are most important for accretion disk structure around Schwarzschild and 
Kerr black holes (see Artemova, Bjornsson & Novikov 1996; Das 2002, and references therein for further discussion). 

Introduction of such potentials allows one to investigate the complicated physical processes taking place in disc 
accretion in a semi- Newtonian framework by avoiding pure general relativistic calculations so that most of the features 
of spacetime around a compact object are retained and some crucial properties of the analogous relativistic solutions of 
disc structure could be reproduced with high accuracy. Hence, those potentials might be designated as 'pseudo-Kerr' 
or 'pseudo- Schwarzschild' potentials, depending on whether they are used to mimic the space time around a rapidly 
rotating or non rotating/ slowly rotating (Kerr parameter a ~ 0) black holes respectively. Below we describe four 
such pseudo Schwarzschild potentials on which we will concentrate in this article. In this section, as well as in the 
following sections, we will use the value of to be equal to IGMbh I <? ■ 

It is important to note that as long as one is not interested in astrophysical processes extremely close (within 
1 — 2 Tg) to a black hole horizon, one may safely use the following black hole potentials to study accretion on to a 
Schwarzschild black hole with the advantage that use of these potentials would simplify calculations by allowing one 
to use some basic features of flat geometry (additivity of energy or de-coupling of various energy components, i.e., 

thermal (:^^), Kinetic (^) or gravitational ($) etc., see subsequent discussions) which is not possible for calculations 

in a purely Schwarzschild or a Kerr metric. Also, one can study more complex many body problems such as accretion 
from an ensemble of companions or overall efficiency of accretion onto an ensemble of black holes in a galaxy or for 
studying numerical hydrodynamic accretion flows around a black hole etc. as simply as can be done in a Newtonian 
framework, but with far better accuracy. So a comparative study of multi-transonic accretion flow using all these 
potentials might be quite useful in understanding some important features of the analogue properties of astrophysical 
accretion. 

Also, one of the main 'charms' of the classical analogue gravity formalism is that even if the governing equations 
for fluid flow is completely non-relativistic (Newtonian), the propagation of acoustic fluctuations embedded into it 
are described by a curved pseudo-Riemannian geometry. In connection to astrophysical accretion, one of the best 
ways to manifest such interesting effect would be to study the analogue effects in the Newtonian and post-Newtonian 
accretion flow. However, one should be careful in using these potentials because none of these potentials discussed 
in the subsequent paragraphs are 'exact' in a sense that they are not directly derivable from the Einstein equations. 
These potentials could only be used to obtain more accurate correction terms over and above the pure Newtonian 
results and any 'radically' new results obtained using these potentials should be cross-checked very carefully with the 
exact general relativistic theory. 

Paczyhski and Wiita (1980) proposed a pseudo-schwarzschild potential of the form 



which accurately reproduces the positions of the marginally stable orbit Vg and the marginally bound orbit r;,, and 
provides the value of efficiency to be —0.0625, which is in close agreement with the value obtained in full general 
relativistic calculations. Also the Keplarian distribution of angular momentum obtained using this potential is exactly 
same as that obtained in pure Schwarzschild geometry. It is worth mentioning here that this potential was first 
introduced to study a thick accretion disc with super Eddington Luminosity. Also, it is interesting to note that 
although it had been thought of in terms of disc accretion, $i is spherically symmetric with a scale shift of Vg. 

To analyze the normal modes of acoiistic oscillations within a thin accretion disc around a compact object (slowly 
rotating black hole or weakly magnetized neutron star), Nowak and Wagoner (1991) approximated some of the 
dominant relativistic effects of the accreting black hole (slowly rotating or non-rotating) via a modified Newtonian 
potential of the form 



<I>2 has correct form of Vg as in the Schwarzschild case but is unable to reproduce the value of r-jy. This potential has the 
correct general relativistic value of the angular velocity fig at Vg . Also it reproduces the radial epicyclic frequency 




(122) 
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(for r > Ts) close to its value obtained from general relativistic calculations, and among all black hole potentials, $2 
provides the best approximation for fig and v^- However, this potential gives the value of efficiency as —0.064 which 
is larger than that produced by $1, hence the disc spectrum computed using $2 would be more luminous compared 
to a disc structure studied using <I>i. 

Considering the fact that the free-fall acceleration plays a very crucial role in Newtonian gravity, Artcmova, 
Bjornsson & Novikov (1996) proposed two different black hole potentials to study disc accretion around a non- 
rotating black hole. The first potential proposed by them produces exactly the same value of the free-fall acceleration 
of a test particle at a given value of r as is obtained for a test particle at rest with respect to the Schwarzschild 
reference frame, and is given by 

<^>3 = -I + (l - (123) 

The second one gives the value of the free fall acceleration that is equal to the value of the covariant component of 
the three dimensional free-fall acceleration vector of a test particle that is at rest in the Schwarzschild reference frame 
and is given by 




(124) 



Efficiencies produced by $3 and <I>4 are —0.081 and —0.078 respectively.The magnitude of efficiency produced by $3 
being maximum, calculation of disc structure using $3 will give the maximum amount of energy dissipation and the 
corresponding spectrum would be the most luminous one. Hereafter we will refer to all these four potentials by $i 
in general, where {i = 1,2,3,4} would correspond to <i>i (I12ip . $2 (|122p . $3 (|123p and $4 (I124p respectively. One 
should notice that while all other $i have singularity at r = r^, only $2 has a singularity at r = 0. 



In figure 8 (reproduced from Das & Sarkar 2001), we plot various $i as a function of the radial distance measured 
from the accreting black hole in units of r^. Also in the same plot, purely Newtonian potential is plotted. If we now 
define a quantity Si to be the 'relative stiffness' of a potential $i as: 

r 

(that is, Si is a measure of the numerical value of any ith potential at a radial distance r), we find that for r > 2rg, 

S2 < Sn < Si < S3 < S4, 

which indicates that while $2 is a 'flatter' potential compared to the pure Newtonian potential <I>im, all other pseudo 
potentials are 'steeper' to <I>n for r > 2rg. 

One can write the modulus of free fall acceleration obtained from all 'pseudo' potentials except for $2 in a compact 
form as 



1 



where Si ^ 2, S3 ^ ^ and S4 = 1. 



denotes the absolute value of the space derivative of <I>i, i.e.. 



(125) 



dr 



whereas acceleration produced by $2 can be computed as 

$2 
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(126) 



For axisymmetric accretion, at any radial distance r measured from the accretor, one can define the effective 
potential ^^^'^ (r) to be the summation of the gravitational potential and the centrifugal potential for matter accreting 
under the influence of ith pseudo potential. ^1'^'^ {r) can be expressed as: 



(127) 
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FIG. 8: Newtonian potential and other pseudo-potentials $i(r) {i — 1, 2, 3, 4) are plotted as a function of the logarithmic radial 
distance from the accreting black hole. This figure is reproduced from Das & Sarkar 2001. 



where A(r) is the non-constant distance dependent specific angular momentum of accreting material. One then easily 
shows that A(r) may have an upper limit: 

Ar(r)=riy$iw (128) 

where represents the derivative of $i(r) with respect to r. For weakly viscous or inviscid flow, angular momentum 
can be taken as a constant parameter (A) and (|127p can be approximated as: 

^t"{r)^Mr) + ^ (129) 

For general relativistic treatment of accretion, the effective potential can not be decoupled in to its gravitational and 
centrifugal components. The general relativistic effective potential ^G^{r) (excluding the rest mass) experienced by 



39 



the fluid accreting on to a Schwarzschild black hole can be expressed as: 



r - 1 
r3 - A2 (1 + r) 



- - Jl^^TTT^^ - 1 (130) 



One can understand that the effective potentials in general relativity cannot be obtained by linearly combining its 
gravitational and rotational contributions because various energies in general relativity are combined together to 
produce non-linearly coupled new terms. 



In Fig. 9. (reproduced from Das 2002), we plot ^^•'{r) (obtained from (|129p ) and ^Qjiir) as a function of r in 
logarithmic scale. The value of A is taken to be 2 in units of 2GM/c. ^'^^^ curves for different <i>i are marked 
exclusively in the figure and the curve marked by represents the variation of ^'q^{t') with r. One can observe 
that (r) is in excellent agreement with $^^(f). Only for a very small value of r (r^rg), starts deviating 

from '^'q^{t) and this deviation keeps increasing as matter approaches closer and closer to the event horizon. All other 
^1^^ {r)s approaches to ^'^^('') at a radial distance (measured from the black hole) considerably larger compared to 
the case for <^\^^ (r). If one defines A^'^^(r) to be the measure of the deviation of '^l^^ {r) with ^'q^^t) at any point 
r, 

Af^(r) = $f^(r)-$-Ji(r) 

One observes that A^^''^(r) is always negative for ^l^-^ {r), but for other ^'^■^■^ {r), it normally remains positive for low 
values of A but may become negative for a very high value of A. If |A^''^^(r)| be the modules or the absolute value 
of Aff{r), one can also see that, although only for a very small range of radial distance very close to the event 
horizon, A^'^'^{r) is maximum, for the whole range of distance scale while $1 is the best approximation of general 
relativistic space time, $2 is the worst approximation and $4 and <I>3 are the second and the third best approximation 
as long as the total effective potential experienced by the accreting fluid is concerned. It can be shown that |A^''^^(7')| 
nonlinearly anti-correlates with A. The reason behind this is understandable. As A decreases, rotational mass as well 
as its coupling term with gravitational mass decreases for general relativistic accretion material while for accretion in 
any ^i, centrifugal force becomes weak and gravity dominates; hence deviation from general relativistic case will be 
more prominent because general relativity is basically a manifestation of strong gravity close to the compact objects. 
From the figure it is clear that for ^Q^ir) as well as for all ^l^'^ (r), a peak appears close to the horizon. The height 
of these peaks may roughly be considered as the measure of the strength of the centrifugal barrier encountered by 
the accreting material for respective cases. The deliberate use of the word 'roughly' instead of 'exactly' is due to the 
fact that here we are dealing with fluid accretion, and unlike particle dynamics, the distance at which the strength 
of the centrifugal barrier is maximum, is located further away from the peak of the effective potential because here 
the total pressure contains the contribution due to fluid or 'ram' pressure also. Naturally the peak height for <&^/(r) 
as well as for {r) increases with increase of A and the location of this barrier moves away from the black hole 
with higher values of angular momentum. If the specific angular momentum of accreting material lies between the 
marginally bound and marginally stable value, an accretion disc is formed. For inviscid or weakly viscous flow, the 
higher will be the value of A, the higher will be the strength of the centrifugal barrier and the more will be the amount 
of radial velocity or the thermal energy that the accreting material must have to begin with so that it can be made to 
accrete on to the black hole. In this connection it is important to observe from the figure that accretion under ^i{r) 
will encounter a centrifugal barrier farthest away from the black hole compared to other $i. For accretion under all 
^iS except $i,the strength of centrifugal barrier at a particular distance will be more compared to its value for full 
general relativistic accretion. 

In subsequent sections, we will use the above mentioned potentials to study the analogue effects in spherically 
symmetric and in axisymmetric black hole accretion. 



XV. NEWTONIAN AND POST-NEWTONIAN SPHERICAL ACCRETION AS AN ANALOGUE 

MODEL 



In this section, we study the analogue gravity phenomena in the spherical accretion onto astrophysical black holes 
under the influence of Newtonian as well as various post-Newtonian pseudo-Schwarzschild potentials described above. 
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FIG. 9: i^Btii''')) for pseudo-general relativistic($^^'' (r)) accretion discs as a function of the distance (measured from 

the event horizon in units or Vg) plotted in logarithmic scale. The specific angular momentum is chosen to be 2 in geometric 
units. The figure is reproduced from Das 2002. 

We use the expressions 'post-Newtonian' and 'pseudo-Schwarzschild' synonymously. Our main goal is to provide a 

self-consistent calculation of the analogue horizon temperature Tah in terms of the minimum number of physical 
accretion parameters, and to study the dependence of Tah on various flow properties. This section is largely based 
on Dasgupta, Bilic & Das 2005. 
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A. Equation of Motion 

The non-relativistic equation of motion for spherically accreting matter in a gravitational potential denoted by $ may 
be written as 

du ^ du ^\dp ^ 9$ ^ (131) 
dt dr p dr dr ' 

The first term in (|13ip is the Eulerian time derivative of the dynamical velocity, the second term is the 'advective' 
term, the third term is the momentum deposition due to the pressure gradient and the last term is the gravitational 
force. Another equation necessary to describe the motion of the fluid is the continuity equation 

To integrate the above set of equations, one also needs the equation of state that specifies the intrinsic properties of 
the fluid. We will study accretion described by either a polytropic or an isothermal equation of state. 



B. Sonic Quantities 

1. Polytropic Accretion 

We employ a polytropic equation of state of the form p ~ Kp'' . The sound speed Cg is defined by 



dp 
dp 



constant entropy 



Assuming stationarity of the flow, we flnd the following conservation equations: 
1) Conservation of energy implies constancy of the specific energy £ 



(133) 



£ 



2 7-1 

2) Conservation of the baryon number implies constancy of the accretion rate M 

M = Anpur'^. 

Equation (I134p is obtained from (|13ip . and (|135p follows directly from (|132p . 

Substituting p in terms of Cg and differentiating (|135p with respect to r, we obtain 

c,(l-7) f^^l 
u r 



(134) 



(135) 



(136) 



where ' denotes the derivative with respect to r. Next we differentiate (I134p and eliminating with the help of (|136p 
we obtain 



2c2/r- $' 
u — c^/u 



One thus finds the critical point conditions as 







2 


r— rc 



(137) 



(138) 



As described in section 12.2, here also the critical point and the sonic points are equivalent, and the location of the 
sonic point is identical to the location of the acoustic horizon due to the assumption of stationarity and spherical 
symmetry. Thus, hereafter we denote rh as the sonic point and the sphere of radius rh as the acoustic horizon. 
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Hereafter, the subscript h indicates that a particular quantity is evaluated at r^. The location of the acoustic horizon 
is obtained by solving the algebraic equation 

£ - \ (^) r.^'n - = 0. (139) 
The derivative u'^ at the corresponding sonic point is obtained by solving the quadratic equation 



(1+7) «)' + 2(7-l) 

+ (27-l)^ + < = 0, (140) 

which follows from ()137|1 in the limit r^rh evaluated with the help of I'Hospital's rule. 

Finally, the gradient of the sound speed at the acoustic horizon is obtained by substituting u'^^ obtained from 
(I140p into equation (|136p at the acoustic horizon 




(141) 




2. Isothermal Accretion 



We employ the isothermal equation of state of the form 

P= P^cip, (142) 

M 

where T is the temperature, R and fi are the universal gas constant and the mean molecular weight, respectively. The 
quantity Cs is the isothermal sound speed defined by 



dp 
dp 



er, (143) 



where the derivative is taken at fixed temperature and the constant Q — kb / {p-ttih) with vtlh — "nip being the mass 
of the hydrogen atom. In our model we assume that the accreting matter is predominantly hydrogen, hence /i ~ 1. 
Now, the specific energy equation takes the form 

£ = y + eTlnp + $, (144) 

whereas the accretion rate is given by (jl35p as before. 

The radial change rate of the dynamical velocity is again given by (|137p . From (|137p and with (|143p we find the 
sonic point condition as 



Uh^^^^^Cs^Vef. (145) 

since Cg does not depend on r. The derivative of u at r/j is obtained from (|137p by making use of I'Hospital's rule as 
before. We find 



where the minus sign in front of the square root indicates accretion (the plus would correspond to a wind solution). 
Note that the quantities in equations (|145p and (|146p are functions of the fluid temperature T only. Hence the isother- 
mal spherical accretion can be essentially described as a one-parameter solution of the hydrodynamical equations, 
parameterized by T. 
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C. Analogue Temperature 

From psp in Newtonian limit, i.e., 



\x^\=9oo^ + (147) 

gives a general expression for the temperature of the analogue Hawking radiation in a spherically accreting fluid in 
the Newtonian as well as in any pseudo-Schwarzschild gravitational potential 



J AH 



2'KKb V 2c^ 



1 



cir 



(148) 



The quantities required to calculate the analogue temperature (|148p are obtained using the formalism presented in 
section 15.2. For polytropic accretion, using equations (|136p - (|14ip one finds 



Tah ^ 2 + <S>h / 2 ^ ^7+1 
V 2 \2-rh<^, 



Kf(^)_(l + ^)$'^ (149) 

y 'h 

where f (7) = (0.000757^ - 5.00157 + 3.00075). The quantities $,1, <i>Jj, and <^'l are obtained by calculating the values 
of various potentials at r?i, and is calculated from p39p for an astrophysically relevant choice of {£,7}. 

Note that if (cj, — u')h is negative, one obtains an acoustic white-hole solution. Hence the condition for the 
existence of the acoustic white hole is 




^)y^f(7)-(l+7)*'^ <0. (150) 

Since 7 and can never be negative, and since and <I>'^' are always real for the preferred domain of {£,7}, 
unlike general relativistic spherical accretion, acoustic white-hole solutions are excluded in the astrophysical accretion 
governed by the Newtonian or post-Newtonian potentials. 

For an isothermal flow, the quantity is zero and using ()146p we find 

(151) 

where r^ should be evaluated using (|145p . Clearly, the fluid temperature T completely determines the analogue 
Hawking temperature. Hence, a spherical isothermally accreting astrophysical black hole provides a simple system 
where analogue gravity can be theoretically studied using only one free parameter. 

For both polytropic and isothermal accretion, for certain range of the parameter space, the analogue Hawking 
temperature Tah may become higher than the actual Hawking temperature Ty, see Dasgupta, Bilic & Das 2005 for 
further details. 

XVI. POST-NEWTONIAN MULTI-TRANSONIC ACCRETION DISC AS ANALOGUE MODEL 

In this section, we will study the analogue gravity phenomena for polytropic (adiabatic) and isothermal rotating, 
advective, multi-transonic accretion disc in various pseudo-Schwarzschild potentials described in section 14. 

A. Flow Dynamics and Accretion Variables at the Critical Point 

1. Polytropic Accretion 

The local half-thickness, hi{r) of the disc for any $i can be obtained by balancing the gravitational force by pressure 
gradient and can be expressed as: 



h,{r) = Cs^Jr/ (7$0 (152) 
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where <i>J = d^i/dr. For a non- viscous flow obeying the polytropic equation of state p = Kp^ , integration of radial 
momentum equation: 

du IdP d 
dr p dr dr 

leads to the following energy conservation equation (on the equatorial plane of the disc) in steady state: 



= 



(153) 



A2 



7-1 2r2 



and the continuity equation: 



— [uprhi{r)] = 
dr 



can be integrated to obtain the baryon number conservation equation: 



M 



The entropy accretion rate :i can be expressed as: 



l\-2 



ucspr^ ($^) 



1 ('^'l 3 



(154) 



(155) 



(156) 



(157) 



One can simultaneously solve (|154p and (|157p for any particular $i and for a particular set of values of {£, A, 7}. For 
a particular value of {£,A, 7}, it is now quite straight-forward to derive the space gradient of the acoustic velocity 



' dc. 



)^ and the dynamical flow velocity (^)j for flow in any particular ith black hole potential <i>i 



dcs 
dr 



7-1 



1 



1 du 



7 + 1/ V2 2r udr 



and, 



3 I £iW 



2c5 



«(7 + l) 



where $/' represents the derivative of $/. Hence the critical point condition comes out to be: 

$-(r) +7$-(r) / A2+r3$^(r) 



1 + 7 



3$;(r)+r$"(r) 



(158) 



(159) 



(160) 



Note that the Mach number at the critical point is not equal to unity, rather: 

Hence, the critical points and the sonic points are not equivalent. One needs to calculate the sonic point, which is the 
location of the acoustic horizon, following the procedure as described in section 13.5. 

For any flxed set of {£, A, 7}, the critical points can be obtained by solving the following polynomial of r: 



E- 



2r2 



27 



3$,(r)+r$, (r) 



= 0. 



(162) 



The dynamical velocity gradient at the critical point can be obtained by solving the following equation for [du/dr)r 



^'^y -2M 



7 + 1 \dr J 



7-1 



3 + 

r $;(r) 



27 



7+1 

6(7-1) f'^i{r)\ 6(27-1) 



du\ 



(i + 7)'V'^^Wy 7(7 + 1)' V'^^W/ 7^(7 + 1)' 



3A2 



(163) 
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Where the subscript (c, i) indicates that the corresponding quantities for any ith potential is being measured at its 
corresponding critical point and = 



2. Isothermal Accretion 



The isothermal sound speed is defined as: 



(164) 



where Q = ^^^^ is a constant, mn^Trip being the mass of the hydrogen atom and kb is Boltzmann's constant. The 



local half-thickness hi{r) of the disc for any 3>i(r) can be obtained by balancing the gravitational force by pressure 
gradient and can be expressed as 



hi{r) = e. 



(165) 



Solution of the radial momentum conservation equation and the continuity equation provides the following two integral 
of motion on the equatorial plane of the isothermal accretion disc 



,,2/ ■\ \2 



Constant 



and 



(166) 



T 

M = Qp{r)u{r)r'^ — 



(167) 



The dynamical flow velocity for a particular value of {£, A} can be expressed as 



du 
dr 



(^0 



2rp$"('-) 



<S>'-(r) 



(168) 



where $" = llp^- Since the flow is isothermal, dcg/dr = everywhere identically. 
The critical point condition can be expressed as: 



eT3 



A2 



\ 2rc 2 {% )^^^^ 



(169) 



Note that the Mach number at the critical point is exactly equal to unity, hence the critical points and the sonic 
points are identical for isothermal accretion disc. Therefore, rc is actually the location of the acoustic event horizon 
Th, and for a specific value of {£,X}, Th can be computed by solving the following equation for r^: 



3 2A2 

+ 



r=rh 

The dynamical velocity gradient at the acoustic horizon can be obtained as: 

2 



(170) 



/ du 



+ 



39^r 3A^ 



(171) 
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B. Multi-transonicity and Shock Formation 



As in the case of general relativistic accretion disc, axisymmetric accretion under the influence of a generahzed 
pseudo-Schwarzschild potential also produces multiple critical/sonic points, both for polytropic as well as for the 
isothermal flow. For polytropic flow, (|162p can be solved to obtain various critical points, and the flow equations can 
be integrated from such critical points to flnd the corresponding sonic points. 

For accretion/ wind solutions under the influence of various ^i, one can define the square of the eigenvalue Q in 
the following way (Chaudhury, Ray & Das 2006): 



4rc$'(rc)|c; 



s |r=rc 



(7 + 1)= 



BC 

(7-l)^-27(l+C) + 27— 



BC 



A|(rc) 



47+(7-l)^ + 27— 



(172) 



where 



A 



= $'(rc) 



3, 6 = 1 



$'"(rc) <i>"(rc) 



$"(re) 



C^A + 3, A^(r) =r^^'{r) 



(173) 



For isothermal flows, a similar expression for the related eigenvalues may likewise be derived. The algebra in this case 
is much simpler and it is an easy exercise to assure oneself that for isothermal flows one simply needs to set 7 = 1 in 
P72p . to arrive at a corresponding relation for Cl^. 

A generic conclusion that can be drawn about the critical points from the form of in (|172|) . is that for a 
conserved pseudo-Schwarzschild axisymmetric flow driven by any potential, the only admissible critical points will be 
saddle points and centre-type points. For a saddle point, > 0, while for a centre-type point, Cl^ < 0. Once the 
behaviour of all the physically relevant critical points has been understood in this way, a complete qualitative picture 
of the flow solutions passing through these points (if they are saddle points), or in the neighbourhood of these points 
(if they are centre- type points), can be constructed, along with an impression of the direction that these solutions can 
have in the phase portrait of the flow, see Chaudhury, Ray & Das (2006) for further detail. 

Application of the above mentioned methodology for finding out the nature of the critical point leads to the 
conclusion that for multi-transonic accretion and wind, the inner critical point r™ and the outer critical point r°"* 
are of saddle type ('X' type), whereas the middle critical point r™*"* is of centre type ('O' type). For mono-transonic 
accretion, the critical point will always be of saddle type and will be located either quite close to the event horizon 
(mono-transonic accretion passing through the 'inner type' critical point) or quite far away from the black hole 
(mono-transonic accretion passing through the 'outer type' critical point). 

Hereafter we will use the notation [Pi] for a set of values of [£, A, 7] for polytropic accretion in any particular 
For all <i>i, one finds a signiflcant region of parameter space spanned by [Vi] which allows the multiplicity of critical 
points for accretion as well as for wind where two real physical inner and outer (with respect to the location of the 
black hole event horizon) saddle type critical points r™ and r°"* encompass one centre type unphysical middle sonic 
point r™''' in between. For a particular $i, if Ai [Pi] denotes the universal set representing the entire parameter space 
covering all values of [Pi], and if Bi [Pi] represents one particular subset of Ai [Pi] which contains only the particular 
values of [Pi] for which the above mentioned three critical points are obtained, then Bi [Pi] can further be decomposed 
into two subsets Ci [Pi] and Vi [Pi] such that: 



[P^] C B, [P,] only for S (r™) > S (r°"*) 
'D^ [P^ C B, [P,] only for S (r^) < S (r°"*) 



(174) 



then for [Pi] G Ci [Pi], we get multi-transonic accretion and for [Pi] G Vi [Pi] one obtains multi-transonic wind. 



For the Paczyhski & Wiita (1980) potential $1, in figure 10 we classify the whole [£, X] parameter space for a fixed 
value of 7 = 4/3. The region marked by I represents the values of [£, X] for which accretion will be mono-transonic 
and will pass through the saddle type inner critical point, whereas the region marked by O represents the values 
of [£, X] for which accretion will be mono-transonic and will pass through the saddle type outer critical point. The 
wedge shaped region bounded by heavy solid lines and marked by A (including the shaded region) represents the 
multi-transonic accretion zone for which {£i, Ai) € [Pi] € Ci [Pi] C Bi [Pi], whereas the wedge shaped region bounded 
by the heavy solid line and marked by W represents the multi-transonic wind and mono-transonic accretion zone for 
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FIG. 10; The complete classification of [£,\] for polytropic accretion in the Paczyriski & Wiita (1980) potential ^i. The value 
of 7 is taken to be equal to 4/3. Mono-transonic regions are marked by I (accretion through the inner sonic point only) and O 
(accretion through the outer sonic point only). The regions marked by A and W represents the multi-transonic accretion and 
wind, respectively. The shaded region represents the collection of [S,X] (for 7 = 4/3) for which the stable Rankine-Hugoniot 
shock solutions are obtained. 



which {Si, Aj) e [Vi] e Vi [Pi] C [Pi]- A similar kind of parameter space division can easily be obtained for other 

$i as well, sec Das 2002 and Chaudhury, Ray & Das 2006 for further detail. 

If shock forms in accretion, then [Pi] responsible for shock formation must be somewhere from the region for 
which [Pi] e Ci [Pi], though not all [Pi] G Ci [Pi] will allow shock transition. One can derive (see Das 2002 for further 
detail) the Rankine-Hugoniot shock condition for the generalized potential $i in the following form which will be 
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satisfied only at the shock location 



1 i-i^i 



(l-7)(^^^j f(fe.+th)-/3i(l + /3i-Pcomp)^' + (l + /3i) ' =0 (175) 

where M is the mass accretion rate as defined in ()156p . £(^ki+th) is the total specific thermal plus mechanical energy 
of the accreting fluid: 



(ki+th) 



Pcomp and O are the density compression and entropy enhancement ratio respectively, defined as Pcomp = iP+/P-) 
and 6 = |^S+/S„^ respectively; Pi — 1 — r^^^'^'> and F = Qpcomp, "+" and "_" refer to the post- and pre-shock 
quantities. The shock strength Si (ratio of the pre- to post-shock Mach number of the flow) can be calculated as: 

Si = Pcomp (1 + /3i) (176) 

Equations (|175p and (jl76p cannot be solved analytically because they are non-linearly coupled. However, one can 
solve the above set of equations using iterative numerical techniques. An efficient numerical code has been developed 
in Das 2002, which takes \Pi\ and ^i as its input and can calculate the shock location r^^ along with any sonic or shock 
quantity as a function of \Pi\ . One obtains a two- fold degeneracy for , and the local stability analysis ensures that 
the shock which forms in between the sonic points r°"* and r™"' is stable for all ^i. Hereafter, we will be interested 
only in such stable shocks and related quantities. 

If [Pi] G J-i [Pi] C d [Pi] represents the region of parameter space for which multi-transonic supersonic flows is 
expected to encounter a Rankine-Hugoniot shock at Vgh, where they become hotter, shock compressed and subsonic 
and will again become supersonic only after passing through ri„ before ultimately crossing the event horizon, then 
one can also define [Pi] G Gi [Pi] which is complement of Ti [Pi] related to Ci [Pi] so that for: 



g^ [Pi_ 



[P^] e [P,] and [P,] i \P,] \ , (177) 



the shock location becomes imaginary in Qi [Pi], hence no stable shock forms in that region. Numerical simulation 
shows that (Molteni, Sponholz & Chakrabarti 1996) the shock keeps oscillating back and forth in this region. One 
anticipates that Gi [Pi] is also an important zone which might be responsible for the Quasi-Periodic Oscillation (QPO) 
of the black hole candidates, and the frequency for such QPO can be computed for all pseudo-Schwarzschild potentials 
(see Das 2003 for further details). 

The wedge shaped shaded region in figure 10 represents the [Pi] e J^i [Pi] C d [Pi] zone, for which steady 
standing stable Rankine-Hugoniot shock forms, while the white region of the multi-transonic accretion (marked by 

A) represents the j^^, [P,] [P,] G C, [P^] and [P,] i T, [p!\^ zone. 

Similarly, solution of (jl70p provides the multi-transonic accretion and wind regions for the isothermal accretion 
in various <I>i. The corresponding shock conditions can also be constructed and can be solved for a particular value 
of [r, A] to find the region of parameter space responsible for the formation for stable shock solutions. See Das, 
Pendharkar & Mitra 2003 for details about the multi-transonicty and shock formation in isothermal accretion disc 
around astrophysical black holes. 



C. Analogue Temperature 

For axisymmetric accretion in Newtonian limit, one obtains (Bilic, Das & Roy 2007) from ([55)1 



1x1- V^=^/(l + 2$)(l-5-2$^ 



(178) 



Hence the analogue temperature for the pseudo-Schwarzschild, axisymmetric, transonic accretion with space depen- 
dent acoustic velocity would be (Bilic, Das & Roy 2007): 



2'nKi 



A2 A2 
(l + 2*)(l--j-2$-j 



1 



(179) 
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As discussed earlier, once the critical points are found by solving (|162p . one can integrate the flow equations to find the 
sonic point r^, which actually is the location of the acoustic horizon r^- One then finds the value of {du/ dr)^^^_^ and 
{dcs/dr)^^^^. Thus once a specific set of values for [£, A, 7] for polytropic accretion is provided, all the corresponding 
terms in (jl79l) could readily be known and one thus comes up with an accurate estimation of Tah, as well as r, the 
ratio of the analogue to the actual Hawking temperature, as a function of [£", A, 7]. 

In figure 11, we demonstrate the variation of t (plotted along the Z axis) on [£, A] (for a fixed value of 7 = 4/3) 
for multi-transonic shocked accretion flow in Paczyhski & Wiita (1980) potential <i>i. [f,A] used to obtain such 
result, corresponds to the shaded region of figure 10 (for which stable Rankine-Hugoniot shock forms in polytropic 
accretion). As discussed in section 13.9, two acoustic black holes are formed at the inner and the outer sonic points, 
and an acoustic white hole is formed at the shock location. The analogue temperature corresponding to the white 
hole in not defined. The red surface in the figure corresponds to the variation of r with [£, A] for the outer acoustic 
horizon (the outer sonic points) and the blue surface corresponds to the variation of r with [£, A] for the inner acoustic 
horizons (the inner sonic points). It is observed that for a fixed value of [£, A, 7], r^in > t^ou*. 

Although the above figure has been obtained for a fixed value of 7 (=4/3), one can obtain the same [r — £ — A] 
variation for any value of 7 producing the multi-transonic shocked accretion flow. In general, r co-relates with 7. 
[r — f — A] variation can also be studied for mono-transonic accretion passing through the inner or the outer sonic 
point only, and for mono-transonic accretion flow in multi-transonic wind region (flow described by [£,A] obtained 
from the W region of the figure 10). 

All the above mentioned variation can also be studied for all other <I>i, see Bilic, Das & Roy (2007) for further 
detail. 



0.12 |- 




FIG. 11: Variation of r on £ and A for multi-transonic shocked accretion in the Paczyhski & Wiita (1980) potential c&i. The red 
surface in the figure corresponds to the variation of r with [£, A] for the outer acoustic horizons and the blue surface corresponds 
to the variation of r with [£, A] for the inner acoustic horizons. This figure has been reproduced from Bilic, Das & Roy (2007) 



It is now easy to calculate the analogue temperature for isothermal axisymmetric accretion in pseudo-Schwarzschild 
potential. Since is a function of the bulk temperature of the flow T, and since for isothermal accretion T is constant 
throughout, the space derivative of the acoustic velocity (dcs/dr) is identically zero everywhere for any potential $i. 
Hence the expression for the analogue temperature can be obtained by setting [dcs/dr) = in (|179p . The dependence 
of Tah on [T, A] has been discussed in Bilic, Das & Roy (2007) in detail. 

XVII. EPILOGUE 

The primary motivation of this review article is to demonstrate the following: 

One can propose that the general relativistic as well as the Newtonian/post- Newtonian accretion flow around an 
astrophysical black hole can be considered as an example of classical analogue gravity model realized in nature. To 
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accomplish this task, one first formulates and solves the equations describing the accre tion processes around black 



holes, and then provides the arguments that such accretion is transonic in general 23l| and the accreting material 
must encounter a sonic point at some specific length scale determined by various accretion parameters. The collection 
of such sonic points forms a null hypersurface, generators of which are the acoustic null geodesies, i.e. the phonon 
trajectories. Such a surface can be shown to be identical with an acoustic event horizon. The acoustic surface 
gravity and the corresponding analogue horizon temperature Tah at the acoustic horizon arc then computed in terms 
of fundamental accretion parameters. Physically, the analogue temperature is associated with the thermal phonon 
radiation analogous to the Hawking radiation of the black-hole horizon. Acoustic white holes can also be generated if 
the accretion flow is multi-transonic and if such multi-transonic black-hole accretion encounters a stable shock. Such 
a white hole, produced at the shock, is always flanked by two acoustic black holes generated at the inner and the 
outer sonic points. 

At this point one might as a crucial question: Does the accretion processes only around a black hole represents an 
analogue system, or any kind of astrophysical accretion exhibits the analogue gravity phenomena in generaU From 
the discussions presented in this article, one understands that two of the essential requirements for a physical system 
to manifest the classical analogue gravity effects are the following: 

1. The system should consists of transonic, barotropic fluid, and the fluid should, preferably, be inviscid in order 
not to violate the Lorentzian invariance. 

2. An acoustic perturbation (or equivalent perturbation, a surface gravity wave for example, see, e.g., Schiitzhold 
& Unruh 2002) should propagate within such fluid for which a space time metric can be constructed. Such 
metric should incorporate a singularity (not always in a formal sense though) , from which one can come up with 
the notion of the acoustic horizon. 

Hence, it is obvious that hydrodynamic, non-dissipative accretion onto any astrophysical object should manifest the 
analogue gravity phenomena, if such accretion exhibits transonic properties, and if such accreting fluid configuration 
possesses a specific well defined symmetry (spherically symmetric or axisymmetric flow, for example). Hence, hydro- 
dynamic, transonic, astrophysical accretion possessing a suitable symmetric geometrical configuration may exhibit 
the analogue properties in general, where the accretor resembles the sink. 

Transonic accretion in astrophysics can he conceived to constitute an extremely important class of classical analogue 
gravity model. Among all the classical analogue systems studied in the literature so far, only an accreting astrophysical 
object incorporates gravity (through the general body force term in the Euler's equation, even if the accretion is studied 
within the framework of the Newtonian space-time) in the analogue model. Also, the simplest possible analogue 
model may be constructed for such objects. For example, the spherically symmetric accretion of isothermal fluid onto 
a Newtonian/semi- Newtonian gravitating mass constitutes an analogue system which can be completely determined 
using a single parameter, the bulk flow temperature of the infalling material (see section 15.3). 

However, among all the accreting astrophysical systems capable of manifesting the classical analogue effect, black 
hole accretion process deserves a very special status. The accreting astrophysical black holes are the only real physical 
candidates for which both the black-hole event horizon and the analogue sonic horizon may co-exist. Hence, the 
application of the analogue Hawking effect to the theory of transonic black hole accretion will be useful to compare 
the properties of these two types of horizons. 

Recently, the spacetime geometry on the equatorial slice through a Kerr black hole has been shown to be equiv- 
alent to the geometry experienced by phonons in a rotating fluid vortex (Visser & Weinfurtner 2005). Since many 
astrophysical black holes are expected to possess non-zero spin (the Kerr parameter a), a clear understanding of the 
influence of spin on analogue models will be of great importance. Some important features on the dependence of the 
analogue temperature on the black hole spin angular momentum of an astrophysical black hole has been discussed in 
this article. In section 13.9 (Fig. 7 and related discussions), it has been shown that the black hole spin does influence 
the analogue gravity effect in a rotating relativistic fluid around it. Also the spin (of the black hole) - angular mo- 
mentum (of the accreting material) coupling modulates such effect. Analogue effect is more prominent for retrograde 
(counter-rotating) flow, resulting a higher value of the corresponding analogue temperature. 

In connection to the acoustic geometry, one can define an 'anti-trapped surface' to be a hypersurface in which 
the fluid flow will be outward directed with the normal component of the three-velocity greater than the local speed 
of sound. In stationary geometry, an anti-trapped surface will naturally be constructed by the collection of sonic 
points corresponding to a spherically symmetric or axisymmetric transonic wind solution emanating out from an 
astrophysical source. Transonic outflow (wind) is ubiquitous in astrophysics, spanning a wide range from solar/stellar 
winds to large-scale outflows from active galaxies, quasars, galactic micro-quasars and energetic gamma ray bursts 
(GRB). In section 13.7, it has been shown how to identify the critical and the sonic points corresponding to the wind 
solutions. Such a scheme can be useful in studying the transonic properties of outflow from astrophysical sources. 
Hence the formalism presented in this paper can be applied to study the analogue effects in transonic winds as well. 
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Recently Kinoshita, Sendouda & Takahashi (2004) performed the causality analysis of the spherical GRB outflow 
using the concept of effective acoustic geometry. Such an investigation can be extended into a more robust form by 
incorporating the kind of work presented in this article, to study the causal structure of the transonic GRB outflows 
in axisymmetry, i.e. for energetic directed outflow originating from a black-hole accretion disc system progenitor. 

In connection to the study of accreting black hole system as a classical analogue gravity model, so far the analogy 
has been applied to dcsc;ribe the classical perturbation of the fluid in terms of a field satisfying the wave equation 
in an effective geometry. Such works do not aim to provide a formulation by which the phonon field generated in 
this system could be quantized. To accomplish this task, one would need to show that the effective action for the 
acoustic perturbation is equivalent to a field theoretical action in curved space, and the corresponding commutation 
and dispersion relations should directly follow (see, e.g., Unruh & Schiitzhold 2003). Such considerations are beyond 
the scope of this article. 

While describing the accretion disc dynamics, the viscous transport of the angular momentum is not explicitly 
taken into account. Viscosity, however, is quite a subtle issue in studying the analogue effects for disc accretion. Even 
thirty three years after the discovery of standard accretion disc theory (Shakura & Sunyaev, 1973; Novikov & Thorne 
1973), exact modeling of viscous transonic black- hole accretion, including proper heating and cooling mechanisms, is 
still quite an arduous task, even for a Newtonian flow, let alone for general relativistic accretion. On the other hand, 
from the analogue model point of view, viscosity is likely to destroy Lorenz invariance, and hence the assumptions 
behind building up an analogue model may not be quite consistent. Nevertheless, extremely large radial velocity close 
to the black hole implies rj„/ <C r^iso where rj„/ and Tyisc are the infall and the viscous time scales, respectively. 
Large radial velocities even at larger distances are due to the fact that the angular momentum content of the accreting 
fluid is relatively low (Beloborodov & lUarionov 1991; Igumenshchev & Beloborodov 1997; Proga & Begelman 2003). 
Hence, the assumption of inviscid flow is not unjustified from an astrophysical point of view. However, one of the 
most significant effects of the introduction of viscosity would be the reduction of the angular momentum. It has been 
observed that the location of the sonic points anti-correlates with A, i.e. weakly rotating flow makes the dynamical 
velocity gradient steeper, which indicates that for viscous flow the acoustic horizons will be pushed further out and 
the flow would become supersonic at a larger distance for the same set of other initial boundary conditions. 

In section 13.2, while constructing the geometry of the general relativistic accretion disc, the expression for the 
disc height has been derived using the prescription of Abramowicz, Lanza & Percival (1997). However, a number 
of other models for the disc height exist in the literature (Novikov & Thorne 1973; Riffert & Herold 1995; Pariev 
1996; Peitz & Appl 1997; Lasota & Abramowicz 1997) The use of any other disc height model would not alter our 
conclusion that black- hole accretion disc solutions form an important class of analogue gravity models (see, e.g.. Das 
2004 for further details about the investigation of the relativistic disc dynamics using the disc height proposed by 
Lasota & Abramowicz (1997)). However, the numerical values of Tah and other related quantities would be different 
for different disc heights. 

For all types of accretion discussed here, the analogue temperature Tah is many orders of magnitude lower 
compared with the fluid temperature of accreting matter. However, the study of analogue effects may be measurably 
significant for accretion onto primordial black holes because the analogue as well as the actual Hawking temperature 
may be considerably high for such situations. There may be a possibility that intense Hawking radiation may not allow 
any accretion due to the domination of strong radiation pressure. However, the situation may be completely different 
for Randall- Sundrum type II cosmology, where during the high energy regime of braneworld cosmology, accretion may 
have significant effects on increasing the mass of the primordial black holes (Guedens, Clancy & Liddle 2002; Guedens, 
Clancy & Liddle 2002a; Majumdar 2003). In braneworld scenario, the accretion onto the primordial black holes from 
surrounding radiation bath may completely dominate over the evaporation process as long as radiation dominations 
persists. It would be interesting to investigate the analogue effects in primordial black hole accretion in Randall- 
Sundrum type - II cosmology, to study whether analogue radiation can really dominate over the accretion phase, 
resulting the enhancement of the black hole evaporation process. One may also like to investigate whether the first 
'black hole explosions' due to Hawking radiation would be acoustic-mediated explosions of the medium surrounding 
the primordial black holes. 

In recent years, considerable attention has been focused on the study of gravitational collapse of massive matter 
clump, in particular, on the investigation of the final fate of such collapse (for a review see, e.g., Krolak 1999). 
Goswami and Joshi (2004) have studied the role of the equation of state and initial data in determining the final fate 
of the continual spherical collapse of barotropic fiuid in terms of naked singularities and the black-hole formation. It 
is tempting to study the analogue effects in such a collapse model. Since at some stage the velocity of the collapsing 
fluid will exceed the velocity of local acoustic perturbation one might encounter a sonic horizons at the radial locations 
of the corresponding transonic points in a stationary configuration. One should, however, be careful about the issue 
that many results in analogue models are based on the assumption of a stationary flow, whereas a collapse scenario 
is a full time dependent dynamical process. 

The correspondence between general relativity and analogue gravity has so far been exploited only on a kinemat- 
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ical, i.e. geometrical level. The analogue gravity systems lack a proper dynamical scheme, such as Einstein's field 
equations in general relativity and hence the analogy is not complete. A certain progress in this direction has recently 
been made by Cadoni and Mignemi (Cadoni 2005; Cadoni & Mignemi 2005), who have established a dynamical 
correspondence between analogue and dilaton gravity in 1+1 dimensions. We believe that the approach presented in 
this article in which an arbitrary background geometry serves as a source for fluid dynamics may shed a new light 
towards a full analogy between general relativity and analogue gravity. 
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